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HIGHLIGHT  SUMI-IARY 

The  phenomena  of  frost  heave  and  loss  of  strength  on  thawing 
produce  serious  highway  damage  in  many  areas  of  the  United  States.  The 
frost  susceptibility  of  soils  is  commonly  predicted  by  grain  size  and 
textural  criteria,  which  can  not  take  into  account  the  compacted  soil 
conditions.  A  parameter  which  does  change  with  the  compaction  variables 
is  pore  size  distribution,  and  the  present  study  focused  on  developing 
experimental  relations  between  the  distribution  of  porosity  and  frost 
heaving  rate. 

For  this  study,  rapid  freezing  tests  were  conducted  to  evaluate 
the  heaving  rate  of  three  different  mixtures  of  an  Indiana  silt  and 
Edgar  plastic  kaolin.  Each  soil  mixture  was  compacted  at  two  or 
three  different  energy  levels,  and  at  three  moisture  levels,  viz., 
wet  of  optimum,  optimum,  and  dry  of  optimum. 

To  obtain  pore  size  distributions  of  the  compacted  soils, 
mercury  intrusion  tests  were  performed.   Since  this  procedure  requires 
the  soil  to  be  free  of  moisture,  freeze  drying  techniques  were  used. 
This  method  almost  eliminated  the  volume  change  and  structural 
modification  expected  from  air  or  oven  drying. 

To  relate  frost  heave  to  pore  size  distribution,  the  method 
of  linear  regression  was  used.  The  prediction  equation  ultimately 
selected  had  an  R  value  of  82^,  and  for  the  soils  tested  read: 


xil 


Y  =  -  5.46  -  29.46 


r    h.O 


^0  "  ^0.4 


+  581.1  (X^  q) 


where  Y    =  heave  rate  in  ram/day, 

X  _  =»  cumulative  porosity  for  pores  >  3.0  ym  but  <  300  ym, 

X^ '   =  total  cumulative  porosity, 

X_  ,  =  cumulative  porosity  for  pores  >  0.4  ym. 

The  above  equation  applies  to  all  test  soil  types  and  tested  levels  of 

compaction  variables. 

In  physical  explanation  of  the  above  equation,  it  is  believed  that 
the  large  pores,  quantified  by  the  value  X»  .,  offer  the  least  resistance 
for  water  to  move  to  the  freezing  front,  whereas  the  smaller  pores 
represent  the  amount  of  available  energy  to  do  the  work  of  frost  heave. 
For  soils  of  the  same  type,  the  fine  pores  are  essentially  independent 
of  compaction  and  water  contents  and  thus  the  larger  pores  will  control 
the  frost  heave.   However,  when  evaluating  different  soils,  heave  will 
be  controlled  by  both  fine  and  large  pores. 


INTRODUCTION 

In  many  areas  of  the  United  States,  the  heaving  of  pavement  due  to 
freezing  soil  is  a  perplexing  problem.   Serious  heaving  of  soils  upon 
freezing  is  not  due  to  the  freezing  of  the  local  pore  water  in  the  soil, 
but  to  the  phenomenon  of  water  migration  to  the  freezing  front  and 
formation  of  ice  lenses.  Loss  of  soil  support  on  thawing  is  a  problem 
which  varies  roughly  with  the  amount  of  heave.   For  significant  frost 
heave  to  occur,  the  soil  must  be  frost  susceptible,  proper  thermal 
conditions  must  exist,  and  a  water  table  reasonably  close  to  the  freezing 
front  must  be  present.   The  major  concern  in  this  study  is  the  first  of 
these  requirements,  that  of  soil  frost  susceptibility,  and  how  its  degree 
may  be  better  predicted  through  correlation  with  the  distribution  of  soil 
porosity. 

Soils  are  directly  evaluated  as  to     frost  susceptibility  by  their 
rates  of  heaving  in  freezing  tests.  More  approximate  predictions  are 
based  upon  the  soil  texture  or  grain  size  distribution.   The  latter 
approach  can  not  take  into  account  variations  in  density  and  soil  fabric. 
Furthermore,  this  method  may  be  overly  conservative  for  clayey  soils. 
Due  to  the  increasing  use  of  clayey  and  potentially  frost  susceptible 
soils  for  highways,  criteria  better  than  the  current  grain  size  ones  are 
needed . 

Pore  size  controls  the  migration  of  water  in  soil  and  the  mechaniSB 
of  frost  heave  to  a  great  degree.  For  this  reason,  it  was  felt  that  a 
study  on  the  relation  of  frost  heave  to  pore  size  distribution  would  be 
fruitful. 


The  determination  of  pore  size  distribution  was  accomplished  with 
the  method  of  mercury  Intrusion,  which  has  been  successfully  used  with 
soils  In  prior  studies  at  Purdue  by  Ahmed,  Lovell  and  Diamond  (1974)  and 
Bhasln  (1975) .   Since  mercury  Intrusion  requires  the  soil  to  be  dry, 
freeze  drying  techniques  were  used  to  minimize  changes  In  the  soil 
structure  during  Its  drying  process. 

As  an  Indication  of  the  ability  of  the  soil  to  heave  during 
freezing,  rapid  freezing  tests  were  performed.  The  test  method  was 
developed  by  Kaplar  (1968)  and  Zoller  (1973),  however,  for  this  study 
modified  procedures  and  apparatus  were  developed  to  better  suit  the 
clayey  silt  test  soils. 

To  determine  the  relation  of  frost  susceptibility  to  pore  size 
distribution,  the  method  of  linear  regression  was  used.   This  method 
allows  a  prediction  equation  to  be  developed,  along  with  statistical 
measures  of  its  reliability. 


LITERATURE  REVIEW 

Frost  Heave  Theory 

The  theoretical  reason  and  mechanism  for  soil  heave  upon  freezing 
haie  been  looked  at  for  many  years.   However,  due  to  the  complexity  of 
the  chemical,  thermal  and  molecular  process  of  moisture  migration  during 
freezing,  a  completely  accurate  mathematical  theory  has  yet  to  be 
developed.   There  are  however,  certain  facets  of  frost  heave  theory 
that  are  widely  accepted,  and  these  will  be  discussed  first. 

The  phenomenon  of  ice  lens  formation  in  frozen  soils  has  been 
observed  for  many  years.   Beskow  (1935)  states  that  Runeberg  observed 
clean  ice  in  frozen  soil  in  1765.  Another  early  work  concerning  frost 
heave  was  the  "capillary  theory  of  freezing  pores"  by  Shtukenberg  In 
1885.  Tsytovich  (1975)  states  that  Shtukenberg  first  recognized  that 
frost  heaving  is  caused  by  migration  of  moisture  toward  the  freezing 
front,  and  proposed  the  above  theory  as  an  explanation. 

An  early  concept  in  the  modem  explanation  of  frost  heave  is  that 
a  film  of  water  separates  the  ice  phase  from  the  soil  grains,  and  that 
the  movement  of  water  to  feed  ice  lenses  is  through  these  thin  films. 
This  was  postulated  by  Taber  (1929)  and  Beskow  (1935)  around  1930,  and 
is  still  fundamental  to  most  water  migration  theories.   Corte  (1963) 
showed  experimentally  that  such  a  film  does  exist  by  freezing  water 
from  the  bottom  upward  in  a  transparent  cylinder.   Once  the  water  was 
partially  frozen,  Corte  sprinkled  several  size  ranges  of  small  soil 
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particles  on  top  of  the  ice.  As  freezing  continued,  some  particles  were 
engulfed  and  others  moved  upward  with  the  advancing  front. 

A  second  concept  fundamental  to  most  frost  heave  theories  is  the 
role  of  freezing  point  depression.  The  fact  that  large  amounts  of 
unfrozen  water  exist  in  colloidal  soil  can  be  shown  by  calorimetric 
studies.  Lovell  (1957a)  found  that  even  at  temperatures  as  low  as  -25  C, 
only  62%  of  the  original  moisture  was  frozen  for  a  silty  clay  compacted 
at  Standard  Proctor  optimum  moisture  content. 

Schofield  (1935)  found  that  as  the  water  content  of  the  soil 
decreased  during  freezing,  the  freezing  point  was  depressed  further  below 
0  C.  He  also  found  that  as  the  water  content  decreased,  the  soil  suction 
increased.  He  therefore  proposed  the  equation. 


H  =  r;^AT  (1) 

T  g 


where  U  =  suction, 

L  =  latent  heat  of  freezing, 
AT  =  change  in  temperature 
g  =  gravitational  acceleration, 

for  the  relation  of  pressure  deficiency  at  a  given  water  content  and  the 

temperature  when  the  water  begins  to  freeze. 

Williams  (1968)  performed  suction  tests  to  give  suction-moisture 
content  relations  and  calorimeter  studies  for  moisture  content  versus 
temperature  depression.  The  results  agreed  closely  with  those  of 
Schofield. 

The  freezing  point  depression  is  often  attributed  to  the  curvature 
of  the  ice-water  interface.   Osier  (1967)  gives  the  formula  in  finite 
difference  form  as 
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where  AT  =•  freezing  point  depression, 

o^y  =  ice-water  interfacial  energy, 

T  =  freezing  point  with  no  curvature  or  interface, 
p-;^  =  density  of  ice, 
L  =  latent  heat  of  fusion, 
T.^^   =  radius  of  curvature  of  the  ice-water  interface. 

Williams  (1968)  points  out  that  the  equation  is  only  good  for  pores 

large  enough  to  have  substantial  amounts  of  free  water,  since  a.     and 

L  vary  in  a  complex  matter  for  the  "bound"  water  close  to  the  particle 

surface.  He  states  that  at  -1.0  C,  the  adsorbed  water  is  about  7%   of 

Vj    ,   where  r.   is  equal  to  4.98  x  10   cm,  and  the  thickness  of  the 
iw        iw     ^ 

adsorbed  water  was  considered  to  be  about  7.0  x  10   cm.  However,  at 
-3.5  C,  the  adsorbed  water  accounts  for  about  2/3  of  the  calculated 
radius  of  the  pore.   Thus  this  formula  is  not  good  for  predicting  large 
freezing  point  depressions,  as  encountered  in  clay  soils.   The  fact  that 
the  ice-water  interface  concept  has  its  limitations  was  also  brought  up 
by  Kolaian  and  Low  (1963).   They  point  out  that  there  must  be  a 
significant  lowering  of  the  freezing  point  of  adsorbed  water  due  to  ita 
interaction  with  the  mineral  particle. 

A  more  comprehensive  freezing  point  relation  was  given  by  Low, 
Anderson  and  Hoekstra  (1968)  who  related  the  freezing  point  depression 
to  the  partial  molar  free  energy  of  the  soil  water.  Their  equation  is 
rather  involved,  however  with  the  use  of  the  computer,  they  were  able 
to  calculate  freezing  point  depression  as  a  function  of  partial  molar 
free  energy  and  the  negative  of  the  relative  partial  molar  heat  content 
of  the  soil  water. 
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In  relating  freezing  point  depression  to  frost  heave  theory,  Scott 
(1969)  states  that  for  coarse,  relatively  uniform  soils,  the  temperature 
depression  will  be  only  a  few  hundredths  of  a  degree  C.  In  such  a  soil, 
the  pore  water  will  freeze  homogeneously  throughout  the  soil  mass,  at 
near  0  C,  and  as  a  result,  little  heave  will  take  place. 

An  entirely  different  situation  arises  with  fine  grained  soils, 
where  several  degrees  C  of  super  cooling  are  required  before  nuclei  form 
and  grow  in  the  soil  pores.   Once  nucleation  occurs  in  the  pore,  growth 
of  the  ice  cluster  proceeds  until  the  ice  almost  fills  the  entire  pore. 
During  this  time,  the  temperature  in  the  pore  either  remains  the  same  or 
rises,  due  to  the  release  of  latent  heat  on  freezing.  At  this  point, 
nucleation  in  the  next  lower  pore  may  result,  or  the  ice  filling  the 
pore  may  propagate  through  the  small  channel  into  the  next  pore.  Either 
event  requires  a  sufficient  amount  of  super  cooling.  An  alternative  is 
for  moisture  to  migrate  to  the  freezing  pore,  with  the  resulting 
enlargement  of  the  pore  eventually  forming  an  ice  lens. 

The  relation  of  the  geometry  of  the  ice-water  interface  to  frost 
heaving  pressures  has  been  used  in  theories  by  Everett  and  Haynes  (1965), 
Penner  (1973),  Williams  (1968)  and  many  others.  Williams  uses  the 
commonly  accepted  formula. 


2^iw 

^Pi-Pw>=-^  (3) 
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where  p.  =  pressure  on  ice, 

p  =  pressure  on  water, 
w 

He  states  that  when. 
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where  r  =  radius  of  connecting  pore  channel, 

the  ice  front  will  pass  through  the  pores.  However  if, 

t  \  ^    iw  (5) 

(Pj  -  P  )  <  

*^i    w     r 
c 

the  ice  front  can  not  pass  through  the  pore  channels  and  ice  lenses  will 
form.  I 

A  different  approach  to  the  motive  force  for  migration  of  water  is 
the  osmotic  theory  proposed  by  Cass  and  Miller  (1959)  and  others.   This 
theory  is  based  on  the  idea  that  due  to  a  concentration  of  ions  in  an 
electric  double  layer  around  soil,  an  osmotic  pressure  potential  exists. 
The  authors  believe  that  an  ice-water  interface  can  be  considered  as  a 
surface  effectively  midway  between  particles,  and  that  an  osmotic 
pressure  difference  exists  between  the  interface  and  the  less 
concentrated  solution  in  underlying  pores.   Due  to  this  difference, 
underlying  pore  v/ater  flows  to,  and  then  joins  the  ice  lens,  with  heave 
resulting.   Since  ice  forms  as  a  pure  phase  without  electrolytics,  the 
concentration  of  ions  below  the  ice  lens  tends  to  increase,  whereas  flow 
from  pore  water  tends  to  dilute  the  concentration. 

The  authors  also  believe  that  the  concentration  of  ions  result  in  a 
depression  of  the  double  layer  temperature,  with  the  water  freezing  at  a 
lower  temperature.   They  state  that  ice  lenses  do  not  form  in  a  coarse 
soil  because  water  in  the  double  layer  is  insignificant  compared  to  the 
free  water  in  the  much  larger  pores. 


8 

With  regards  to  the  validity  of  the  osmotic  theory,  Tstovich  (1975) 
says  that  Bozhenova  found  that  perceptible  results  of  osmosis  are 
observed  only  at  rather  high  salt  concentrations  in  the  pore  water. 

Another  driving  force  in  the  migration  of  water  to  the  freezing 
front  is  that  of  chemical  action.  Tstovich  (1975)  describes  work  done 
under  his  direction  by  Nersisova,  which  considered  the  influence  of 
physio chemical  properties  of  soil  surfaces  on  the  migration  of  moisture 
and  ice  segregation.   For  the  study,  kaolin  and  montmorillonite  clays 
were  saturated  with  various  exchange  cations.   It  was  found  that  greater 
migration  of  water  and  heaving  were  observed  when  the  clays  were 
saturated  with  multivalent  cations  and  less  with  univalent  cations. 

An  entirely  different  kind  of  moisture  migration  during  freezing  is 
from  vapor  transportation.   Tsytovich  (1975)  says  that  Lebedev  found 
experimentally  that  migration  of  moisture  from  moist,  but  unsaturated, 
soils  is  due  to  movement  of  vapor  from  locations  of  higher  to  those  of 
lower  vapor  pressure.  However,  Jumikis  (1957)  found  the  amounts  of  soil 
moisture  transferred  upward  to  the  freezing  front  in  the  form  of  vapor 
to  be  insignificant. 

The  basic  adsorbed  water  film  theory  is  probably  the  most  widely 
accepted  and  is  described  in  various  ways  by  Kaplar  (1970),  Hoekstra, 
Chamberlain  and  Frate  (1965),  and  by  Keinonen  (1973).  According  to  this 
theory,  water  migrates  from  thick  films  with  more  mobile  molecules  toward 
thin  films  in  the  partially  ice  filled  pores.   The  reason  for  this  is  that 
Gibbs'  free  energy  is  smaller  for  the  thinner  water  films  and  water  flows 
in  the  direction  of  decreasing  free  energy. 

Keinonen  (1973)  expresses  the  Gibbs'  free  energy  as. 


9r 
G  -  u  +  pv  -  Ts  (6) 

where  G  =  Glbbs'  free  energy, 
u  =  internal  energy » 
p  =  pressure, 
V  =  volume, 
T  =  Temperature, 
8  =  entropy. 

He  states  that  when  the  soil  water  temperature  decreases  below  0  C,  the 

outermost  part  of  the  adsorbed  water  freezes  and  the  f^lm  becomes 

thinner.  When  the  water  molecule  freezes,  it  loses  its  dipolar 

attraction  to  the  soil  particle.  As  a  result,  the  inner  pressure  in  the 

remaining  water  film  decreases,  and  the  Gibbs'  free  energy  consequently 

decreases. 

With  regards  to  actual  frost  heaving,  the  difference  between  the 
free  energy  in  free  ground  water  and  in  freezing  soil  limits  the  total 
mechanical  work  done  during  the  freezing  process.   The  work  in  the 
freezing  process  is  used  both  to  draw  up  water  to  the  freezing  front, 
and  for  frost  heave  pressures. 

Kaplar  (1970)  prefers  to  let  the  soil  physicists  and  physical 
chemists  debate  the  actual  motive  force  for  the  movement  of  water  in 
films.   Instead,  he  looks  at  the  migrating-heaving  system  as  a 
"jacking"  process.   The  thin  films  serve  as  the  jacks,  lifting  the 
overburden  for  the  ice  cluster  to  increase  in  size,  and  in  the  mean 
time  "sucking"  water  to  the  cluster. 

How  involved  one  should  get  in  various  theories  of  frost  heave  and 
water  migration  is  open  to  question,  especially  since  there  appears  to 
be  several  motive  forces  behind  the  process.  Possibly  a  quite  adequate 
explanation  for  the  principle  of  moisture  migration  in  freezing  soils 
Is  given  by  Tsytovich  (1975)  who  says: 
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"...the  migration  of  water  in  freezing  moist  soils  is 
a  process  of  moisture  transport  that  operates 
consistently  whenever  the  equilibrium  state  of  the 
soil's  phases  is  disturbed  and  external  factors 
change. (the  presence  of  temperature,  moisture  content, 
pressure,  mineral  particle  surface  energy,  mobility 
of  molecules  in  water  films,  gradients,  etc)." 

In  summarizing  the  theoretical  and  experimental  work  that  has  been 

accomplished  so  far,  there  is  still  conflict  and  lack  of  understanding 

in  certain  areas,  however,  certain  concepts  are  considered  essentially 

true.   Some  of  these  are: 

1)  Heave  in  fine,  frost  susceptible  soils  is  due  to 
migration  of  moisture  and  not  just  the  slight 
expansion  of  existing  moisture  when  it  freezes. 

2)  Coarse  soils  may  not  heave  appreciably,  due  to 
migrating  forces  being  too  small  to  transport 
v;ater  to  the  freezing  front.   In  some  cases, 
squeezing  of  water  from  the  freezing  front  may 
result. 

3)  Very  plastic  clays  may  only  heave  slightly,  due  to 
excessive  amounts  of  energy  required  to  move  water 
to  the  freezing  front. 

4)  A  film  of  liquid  separates  the  frozen  water  from 
the  soil  particles,  with  the  migration  of  water 

to  the  ice  front  most  likely  being  in  these  films. 

5)  Nucleation  of  water  in  pores  is  partly  a  function 
of  their  curvature  and  size.   In  large  pores, 
nucleation  will  result  at  only  a  few  hundredths  of 
a  degree  C  below  0°C,  whereas  it  may  take  several 
degrees  C  below  0°C  for  nucleation  in  small  pores. 

6)  The  motive  force  for  migration  is  a  complicated 
conglomerate,  which  may  include  chemical,  osmotic, 
vapor  pressure,  and  other  effects. 

7)  Although  the  heaving  mechanism  is  complicated,  it 
is  at  least  partially  dependent  on  the  size  of 
pores,  the  shape  and  size  of  pore  cliannels,  and 
surface  area  of  the  particles. 


11 


Frost  Heave  Tests 

An  important  tool  for  predicting  frost  heave  has  been  the  various 
laboratory  frost  heave  tests.  This  area  is  again  one  in  which  there 
exists  no  complete  agreement,  and  there  seems  to  be  no  universal, 
standard  test  for  all  situations. 

The  U.  S.  Army  CRREL  heave  test  developed  in  the  early  1950' s  is 
probably  the  most  widely  known  one  in  the  U.  S.  As  described  by  Kaplar 
(1974) ,  the  test  is  based  on  one-dimensional  heave  of  samples  frozen 
from  the  top  down  at  a  rate  of  penetration  of  the  0°C  isotherm  of  about 
0.25  to  0.5  inches  per  day.   The  soil  samples  are  frozen  in  acrylic 
tubes,  tapered  so  that  the  inside  diameter  is  5.75  in.  at  the  top  and 
5.5  in.  at  the  bottom.   Before  being  placed  in  the  acrylic  tubes,  the 
soil  is  compacted  in  steel  cylinders  using  the  Providence  Vibrated 
Density  method  for  coarse  soils,  and  AASHTO  test  procedure  T-180-570, 
method  "A"  or  "D",  for  the  fine  grained  soils.   Saturation  is  accomplished 
using  deaired  water  and  a  vacuum. 

A  disadvantage  with  this  method  is  it  is  very  time  consuming,  with 
a  freezing  test  time  of  15  to  20  days.  Also,  the  soil  may  freeze  to  the 
sides  of  the  acrylic  containers,  thus  limiting  heave.   Kapler  (1968) 
believes  that  with  adfreezing  to  the  sides,  a  more  rapid  penetration  rate 
will  result,  due  to  less  latent  heat  release  at  the  freezing  front. 
Since  the  operator  controls  the  rate  of  penetration  of  the  0  C  isotherm, 
once  an  increasing  rate  of  penetration  was  noticed,  he  would  raise  the 
freezing  temperature,  thus  further  decreasing  the  heat  extraction  and 
freezing  rate. 

Based  on  tests  carried  out  by  CRREL,  Casagrande  proposed  the  test 
criteria  shown  in  Table  1. 
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Table  1.   CRREL  Frost  Heave  Classification 

Average  Rate  of  Heave  Relative  Frost  Susceptibility 

(mm  /  day)  Classification 

0-0.5  Negligible 

0.5  -  1.0  Very  Low 

1.0  -  2.0  Low 

2.0  -  4.0  Medium 

4.0  -  8.0  High 

>  8.0  Very  High 

Table  2.   TRRL  Frost  Heave  Classification 

Total  Heave  Relative  Frost  Susceptibility 

(era  for  250  hours)  Classification 

<  1.27  Satisfactory 

1.27  -  1.78  Marginal 

>  1,78  Unsatisfactory 

Table  3.   UNH  Frost  Heave  Classification 

Average  Rate  of  Heave  Relative  Frost  Susceptibility 

(mm  /  day)  Classification 

0-6.5  Negligible 

6.5  -  8.0  Very  Low 

8.0-10.3  Low 

10.3  -  13.0  Medium 

13.0  -  15.0  High 

>  15.0  Very  High 
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The  British  Transportation  and  Road  Research  Laboratory  (TRRL) 
frost  heave  test,  as  described  by  Jacobs  (1965),  is  another  test  widely 
used  for  evaluation  of  freezing  conditions.  As  in  the  CRREL  test, 
freezing  is  unidirectional,  however,  it  is  at  a  constant  tenq)erature. 

The  soil  samples  used  in  the  TRRL  test  are  4  in.  diameter  by  6  in. 
high  cylinders  which  are  i^apped  in  polythene.  The  top  of  the  samples 
are  exposed  to  a  temperature  of  -17  C  and  the  bottoms  rest  on  porous 
ceramic  discs  in  contact  with  water  at  4  C,  with  the  space  between 
samples  filled  with  sand.   The  actual  test  is  carried  out  for  a  period  of 
ten  days,  although  Sutherland  and  Gaskin  (1970)  have  found  experimentally 
that  the  test  time  could  probably  be  reduced  to  about  50  hours. 

The  heave  criterion  for  the  TRRL  test  was  based  on  field  experience 
gained  during  severe  frosts  of  1940  and  1947.  The  frost  susceptibility 
criteria  are  listed  in  Table  2.  ' 

To  compare  the  CRREL  and  TRRL  criteria  and  results,  Sutherland  and 
Gaskin  (1973)  tested  four  mixtures  of  stabilized  pulverized  fly  ash, 
using  both  test  procedures.   Both  tests  classified  the  four  mixes  in  the 
same  increasing  order  of  frost  susceptibility.  Three  of  the  mixes  were 
classified  identically,  with  the  remaining  mix  classified  as  "unsatis- 
factory" by  TRRL  and  "very  low"  by  CRREL. 

Due  to  the  limitations  of  the  CRREL  test,  work  on  a  rapid  freeze 
test  was  performed  by  Kaplar  (1968,  1971)  and  later  Zoller  (1973). 

Kaplar  found  that  with  a  faster  freezing  rate,  heave  rate  increases 
to  a  maximum  and  then  drops  off,  eventually  reaching  a  point  where  heave 
is  due  solely  to  volume  expansion  during  phase  change.  He  therefore 
suggested  that  the  0.25  in.  to  0.5  in.  freezing  penetration  rate  of  the 
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CRUEL  test  was  unnecessarily  slow,  and  more  useful  Information 
concerning  maximum  heave  could  be  obtained  by  using  a  higher  rate  of 
frost  penetration,  along  with  a  constant  temperature. 

Kaplar  also  believed  a  major  problem  with  the  existing  CRREL  test 
was  that  soil  froze  to  the  acrylic  tube,  thus  impeding  heave.  To  reduce 
this  side-wall  friction,  Kaplar  used  a  series  of  acrylic  rings,  designed 
to  separate  during  freezing. 

The  University  of  New  Hampshire  (UNH)  test  developed  by  Zoller 
(1973)  and  his  students  also  is  a  rapid  freeze  test  using  Lucite  rings 
to  reduce  sidewall  friction.  The  test  procedure  used  at  UNH  is  slightly 
different  than  Kaplar 's,  since  a  Peltier  battery  serves  as  a  heat  sink 
during  freezing.   The  battery  is  adjusted  to  a  heat  pumping  capacity 
of  about  65  BTU  per  hour,  with  the  top  of  the  samples  stabilizing  to  a 
temperature  of  about  -4  C. 

A  problem  with  the  use  of  rings  is  that  it  is  difficult  to  saturate 
samples  with  low  permeability.   If  a  vacuum  or  large  gradient  is  applied 
during  saturation,  the  sample  containers  will  leak.   Zoller's  method  of 
saturating  the  sample  by  raising  the  water  level  to  the  top  of  the  sample 
may  work  for  fine  sand  and  some  silts.  However,  it  is  questionable 
whether  a  very  high  degree  of  saturation  can  be  obtained  with  soils 
containing  much  clay. 

The  purpose  of  the  UNH  test  program  was  not  so  much  to  consider  a 
different  laboratory  freezing  method  as  to  correlate  it  with  freezing 
in  the  field.   From  their  testing  program,  the  UNH  criteria  shown  in 
Table  3  were  developed.   Zoller  reports  that  of  the  many  base  course 
materials  tested  during  the  project  life,  no  instances  have  been 
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reported  in  which  materials  classified  as  "low  to  negligible"  in  the 
UNH  rapid  freeze  test  have  subsequently  heaved  detrimentally  in  the 
field. 

Penner  (1972)  has  also  carried  out  experiments  using  a  rapid  freeze 
setup  similar  to  the  UNH  test.  However  Penner  suggests  using  rates  of 
freezing  that  are  related  to  actual  field  conditions,  rather  than  at  one 
arbitrary  applied  temperature. 

In  other  frost  action  tests  carried  out  by  Hoekstra,  Chamberlain, 
and  Frate  (1965),  Osier  (1967),  and  Martin  and  Wissa  (1973),  the  soil 
samples  were  completely  constrained  and  heaving  pressures  were 
measured.  Hoekstra  et  al.  and  others  state  that  the  one-dimensional 
heaving  pressure,  p,  is  given  by 
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where  Oj     =  surface  tension  between  ice  and  water, 

r.  =  effective  radius  of  ice-water  interface, 
iw 

The  researchers  believe  that  heaving  pressures  may  be  a  way  of  predicting 
frost  susceptibility.   However,  the  role  of  permeability  in  transporting 
water  to  the  freezing  front  does  not  seem  to  enter  into  the  heaving 
pressure.  As  a  result,  soils  with  very  low  permeabilities  may  have  high 
heaving  pressures,  yet  heave  very  little. 

Frost  Heave  Criteria 
Before  closing  the  discussion  on  frost  heave,  it  is  worthwhile  to 
discuss  some  of  the  frost  susceptibility  criteria  and  related  factors 
that  are  based  on  soil  parameters. 
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In  the  U.  S.,  probably  the  most  widely  known  frost  susceptibility 
criterion  is  based  on  Casagrande's  (1932)  famous  discussion  in  1932, 
where  he  gave  the  following  rule  of  thumb  for  identifying  frost-sus- 
ceptible soils: 

"Under  natural  freezing  conditions  and  with  sufficient 
water  supply  one  should  expect  considerable  ice  segre- 
gation in  non-uniform  soils  containing  more  than  3  percent 
of  grains  smaller  than  0.02  mm,  and  in  very  uniform  soils 
containing  more  than  10  percent  smaller  than  0.02  mm. 
No  ice  segregation  was  observed  in  soils  containing  less 
than  1  percent  of  grains  smaller  than  0.02  mm,  even  if 
the  ground  water  level  was  as  high  as  the  frost  level." 

A  variation  of  Casagrande's  frost  susceptibility  criterion  is  the 
Corps  of  Engineers  frost  design  classification  system  {Johnson  (1975)}, 
which  is  listed  in  Table  4.  This  test  was  developed  to  account  for  the 
reduced  stability  of  the  various  types  of  frost  susceptible  soils 
during  the  thaw-weakened  period.  Frost  susceptible  soils  are  classified 
into  four  groups  (Fl,  F2,  F3,  F4)  of  increasing  susceptibility  to  frost 
heaving  and/or  weakening  as  a  result  of  frost  melting.  With  the  groups, 
Fl  is  the  least  susceptible,  F4  the  most,  and  F2  and  F3  falling  between 
the  two.   Groups  Fl  and  F2  may  both  experience  equal  ice  segregation, 
however,  Fl  material  may  be  expected  to  show  higher  bearing  capacity 
than  F2  material  during  thaw. 

A  logical  extension  of  textural  and  grain  size  criteria  is  the  pore 
size  distribution  of  the  soil.  Hoekstra,  Chamberlain  and  Frate  (1965) 
stated  that  pore  size  distribution  could  be  a  more  fundamental  parameter. 
They  point  out  that  although  grain  size  and  pore  size  are  related,  this 
relationship  is  obscured  by  factors  such  as  gradation  and  grain  shape. 
As  a  measure  of  pore  size,  the  authors  used  soil  tension  and  were  able 
to  identify  trends  between  the  frost  heaving  process  and  soil  tension. 
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Table  4.   U.S.  Army  Corps  of  Engineers  Frost  Design  Soil  Classification 


Frost 
Group 


Fl 
F2 


Soil  Type 


F3 


F4 


Gravelly  soils 

(a)  Gravelly  soils 

(b)  Sands 

(a)  Gravelly  soils 

(b)  Sands,  except  very 
fine  silty  sands 

(c)  Clays,  PL  >  12 

(a)  All  silts 

(b)  Very  fine  silty  sands 

(c)  Clays,  PL  <  12 

(d)  Varved  clays  and 
other  fine-grained 
banded  sediments 


Percentage 
Finer  Than 
0.02  iran, 
by  weight 


Typical  Soil  Types 
Under  Unified  Soil 
Classification 


3  to  10    GW,  GP,  GW-GM,  GP-GM 


10  to   20 

GM,   GW-GM,   GP-(M 

3  to  15 

SW,    SP,    SM,    SW-SM, 

SP-SM 

>   20 

GM,   GC 

>  15 

SM,    SC 

CL,  CH 

ML,  MH 
>  15    SM 

CL,  CL-CM 

CL  and  ML; 

CL,  ML,  and  SM; 

CL,  CH,  and  ML; 

CL,  CH,  ML,  and  SM 


18 


Csathy  and  Townsend  (1962)  used  a  capillary  rise  test  as  a  basis 
for  determination  of  pore  size  distribution,  and  compared  results  with 
actual  field  frost  behavior.   From  the  comparison,  they  proposed  the 
following  criterion, 

1)  soils  with  P  <  6  are  frost  susceptible,  and 

2)  soils  with  P  >  6  are  non  frost  susceptible, 

where  P^  -  P_-/P^q, 

Pq_  =  the  pore  diameter  such  that  90%  of  the  pores  are  smaller 

than  PgQ. 
P^Q  =  the  pore  diameter  such  that  70%  of  the  pores  are  smaller 

than  P^Q. 

Csathy  and  Townsend  concluded  that  the  pore  size  criterion  was 
considerably  more  efficient  than  the  currently  used  particle  size 
criterion.  However,  Gaskin  and  Raymond  (1973)  point  out  that  Csathy  and 
Townsend  extrapolated  some  of  their  curves  to  obtain  a  P-^  value.   They 
believe  that  this  procedure  is  of  doubtful  validity,  and  state  that  If 
only  curves  that  actually  contain  both  Pq_  and  P^„  values  are  used, 
Csathy  and  Townsend 's  criterion  has  no  real  advantage  over  currently 
used  particle  criterion. 

Gaskin  and  Raymond  (1973)  have  also  looked  at  pore  size  distribution 
from  samples  of  actual  field  soils  and  tried  to  correlate  them  to  the 
field  frost  performance.  However,  instead  of  using  capillary  rise 
methods,  which  take  four  weeks  per  test,  they  used  the  pressure  plate 
suction  test,  along  with  the  mercury  intrusion  test.  Based  on  their 
tests,  the  authors  found  no  correlations  which  were  more  efficient  than 
the  currently  used  criterion  based  on  particle  size. 

Two  objections  might  be  raised  against  Gaskin  and  Raymond's 
experimentation.  The  first  is  the  question  of  disturbance  of  the  pore 
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structure  due  to  sampling  and  drying  that  took  place  before  measurement  of 
the  structure.  The  second  is  that  the  pore  size  distribution  was  only 
carried  down  to  about  the  10  micron  size  range.   For  some  of  the  soils 
tested,  less  than  60%  of  the  pores  were  larger  than  this  size. 
Therefore,  the  actual  complete  pore  size  distribution  was  not  always 
truly  represented. 

When  considering  the  effects  of  density  and  compaction  on  frost 
heave,  Zoller  (1973)  postulates  that  there  is  approximately  a  linear 
relationship  between  heave  rate  and  percent  compaction,  and  that  the 
zero  percent  compaction  intercept  is  approximately  the  same  heave  rate 
value  for  all  materials  tested.  Thus  he  says  that  the  calculated  heave 
rate  at  a  given  percent  compaction  for  any  material  is  a  simple  function 
of  its  heave  rate  at  any  other  compaction  as  expressed  by  the  equation: 

HR  =>  S(PC)  -  15.40  (8) 

where  HR  =  heave  rate,  in  mm, 
PC  =  percent  compaction, 
S  =  slope  of  percent  compaction  -  heave  rate  line. 

The  above  relation  is  somewhat  contrary  to  results  of  an  early  study 

by  Winn  and  Rutledge  (19A0) .  Their  data  indicated  that  for  the  soils 

tested,  there  is  one  density  at  which  frost  heave  occurs  most  readily, 

while  at  higher  or  lower  densities  the  action  is  not  so  pronounced. 

Pore  Size  Determination 

The  present  study  is  based  on  past  work  with  pore  size  done  at 

Purdue  University  and  is  geared  more  toward  the  useful  applications  of 

pore  size  distribution  rather  than  refinement  of  the  technique.  For  this 

reason,  the  present  literature  review  on  pore  size  distribution  is  brief. 
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and  the  reader  is  referred  to  Ahmed  (1971)  and  Bhasin  (1975)  for  a  much 
more  comprehensive  review. 

The  method  used  for  the  present  study  is  mercury  intrusion,  which 
was  first  proposed  by  Washburn  (1921).  Washburn  suggested  that  the 
absolute  pressure,  p,  required  to  intrude  a  cylinderical  pore  of 
diameter,  d,  is  related  by 

4  a  cos  6 

m  , 

P  =  - 1 —  (9) 

where  a  =  surface  tension  of  mercury, 

9  =  contact  angle  between  the  mercury  and  the  pore  wall. 

Since  the  proposal  of  the  above  equation,  there  has  been  much 
discussion  about  the  "correct"  contact  angles  and  surface  tension  values, 
along  with  development  of  practical  mercury  intrusion  apparatus.  However, 
the  basic  equation,  along  with  its  limitations  is  still  used. 

A  primary  limitation  with  mercury  intrusion  is  that  it  measures  the 
diameter  of  the  channel  leading  into  the  pore,  and  not  the  actual  pore 
diameter.  Many  researchers  {Ahmed  (1971),  Orr  (1969)}  point  out  the 
"ink  bottle"  effect,  where  mercury  does  not  flow  into  a  large  pore  with  a 
small  opening  until  the  pressure  is  sufficient  to  intrude  the  pore 
opening.  When  related  to  frost  heave,  this  problem  may  not  be  so  major, 
since  the  advancing  ice  front  and  moisture  migration  to  the  freezing 
front  must  both  pass  through,  and  be  restricted  by,  the  pore  channels. 

A  second  problem  with  the  Washburn  equation  is  that  cylindrical 
pores  are  assumed.  As  pointed  out  by  Bhasin  (1975),  this  assimiption  may 
be  fairly  reasonable  since  differences  due  to  pore  shapes  are  within  an 
order  of  magnitude  while  the  pore  size  range  of  compacted  soils  normally 
extends  over  five  orders  of  magnitude. 
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Sample  Drying  Methods 

For  mercury  intrusion  to  be  effective,  the  soil  sample  must  be 
almost  completely  free  of  moisture.   Since  air  or  oven  drying  often 
produces  excessive  shrinkage,  critical  region  drying  and  freeze  drying 
methods  have  been  developed  to  allow  dehydration  at  essentially  constant 
volume. 

The  critical  region  drying  method  is  based  on  the  idea  that  a  fluid 
at  temperatures  and  pressures  below  the  critical  point  can  occur  in  two 
coexisting  phases,  with  different  specific  volumes.  However,  at 
temperatures  above  the  critical  point,  a  fluid  can  exist  only  as  a  single 
phase  with  no  latent  heat  of  vaporization  and  no  change  in  specific 
volume.   Using  this  concept,  the  critical  region  drying  process  begins  by 
raising  the  temperature  and  pressure  such  that  changes  in  the  specific 
volume  are  slow  in  the  compressed  liquid  region.   Once  the  critical  point 
is  reached,  the  soil  water  is  transformed  into  vapor  and  then  is  removed 
in  the  superheated  vapor  region. 

Using  the  critical  region  drying  method,  Bhasin  (1975)  successfully 
dried  Edgar  plastic  kaolin,  Grundite,  Crosby  silty  clay,  and  a  reddish- 
brown  limestone  residual  clay  with  little  or  no  volume  change.  He  also 
tried  Volclay  bentonite,  but  found  that  the  essentially  montmorillonitic 
clay  shrank  appreciably. 

A  second  method  of  drying  is  the  freeze  drying  method.  With  this 
method,  the  soil  sample  is  quickly  frozen  in  liquid  nitrogen  and  then 
placed  in  a  vacuum  where  the  moisture  in  the  form  of  ice  is  removed  by 
sublimation.  Drying  takes  place  when  the  partial  pressure  of  water 
vapor  at  the  frozen  surface  exceeds  that  of  the  surroundings.  As  Ahmed 
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(1971)  points  out,  the  vacuum  does  not  "suck"  water  out,  but  merely 
reduces  the  concentration  of  gas  molecules  present  in  the  pores  and  thus 
reduces  the  resistance  to  water  vapor  flow. 

Using  the  freeze  drying  process,  Ahmed  dried  Grundite  at  various 
compacted  water  contents.  He  found  that  freeze  drying  produced  less  than 
5%   reduction  of  volume. 

Zimmie  and  Almaleh  (1976)  also  used  freeze  drying,  and  found  from 
some  100  tests  that  the  mean  volume  shrinkage  was  less  than  about  7%   for 
both  a  kaolinite  and  sodium  montmorillonite  clay.   The  clays  were 
compacted  with  a  static  pressure,  at  water  contents  between  10  and  U0% 
for  the  kaolinite  and  120%  to  slightly  over  250%  for  the' montmorillonite. 

Pore  Size  Distribution  of  Compacted  Clays 
From  work  done  by  Sridharan,  Altschaeffl  and  Diamond  (1971),  Ahmed, 
Lovell,  and  Diamond  (1974),  and  Bhasin  (1975),  the  following  conclusions 
can  be  drawn  about  pore  size  distribution  and  its  relation  to  compaction 
variables. 

1)  It  is  generally  agreed  that  most  soils  compacted  on 
the  dry  side  of  optimum  exhibit  greater  space  in 
relatively  large  pores  than  that  present  in  soils 
compacted  at  optimum  or  on  the  wet  side.   Bhasin 
believes  that  this  results  from  dry  side  specimens 
retaining  after  compaction  much  of  the  structure  of 
the  individual  aggregations,  with  large  spaces 
present  between  the  aggregations.  However,  the 
interaggregation  spaces  are  completely  eliminated 

by  the  time  the  water  content  is  increased  to  the     ' 
wet  side. 

2)  Bhasin  found  that  increasing  the  compactive  effort 
on  the  dry  side  of  optimum  water  content  diminished 
the  quantity  of  the  larger  pores,  whereas  increasing 
the  compactive  effort  had  little  effect  on  the  pore 
size  distribution  for  soils  compacted  on  the  wet  side. 
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3)  For  the  Grundlte  samples  tested  by  Ahmed,  the  pore  size 
distributions  were  affected  very  little  when  samples 
were  compacted  to  the  same  moisture-unit  weight  conditions 
by  different  methods  of  compaction. 
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SOILS  STUDIED 

The   soils  were  selected  'basecL  on  projected  frost  heave   and   fac- 
ility of  pore   size   distribution  tests.     With  regards  to   frost  heave, 
silt  generally  is   the  highest  heaving  of  any  soil.     Therefore,   it  was 
decided  to  use   a  relatively  pure  silt   and  observe  the  effects   of 
different  water  contents   and  clay   contents   on   frost  heave   and  pore  size 
distribution.      A  second  factor  favoring  the   choice  of  a  clay-silt   soil 
was   the   requirement  that  the  porosimeter  sample  be  very  small.     Thus, 
soils   containing  particles   larger  than   fine  sand  are  not   really  practical 
with   this    equipment . 

The   silt  was   obtained  from  a  natural   loess   deposit   at  the  second 
bench  level  on  the   east   side  of  US   ^1,   just   south  of  Patoka,    Indiana. 
For  the   frost  heave  and  pore  size  tests   the  silt  was   air  dried  and 
passed  through   a  number  ^0   sieve.      The  silt  was    found  to  have   a  specific 
gravity   of  2. 73,   liquid  limit   of  2k  and  no  plasticity.     The  AASHTO 
cleissification  of  this   type  of  soil   is  A-U.      A  grain  size   distribution 
for   the  natural  silt   is   shown  in  Figure   1. 

The   clay   that  was   mixed  with   the  silt   so  as   to  give  different  soil 
combinations  was   Edgar  plastic  kaolin.     This    clay   is   commercially 
processed  to  remove  material   coarser  than   kO  micrometers,    dried  in  a 
tunnel   dryer  at  temperatures   ranging  from  300°   to  1+50°F,   and  then 
pulverized.      The  manufacturer  states   that  the   clay  mineral   composition 
is   99-5!^  kaolinite,  with  the   remainder  as  micro  sized  alpha  quartz   and 
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mica.  The  specific  gravity  of  the  clay  is  2.65. 

Three  silt-clay  mixtures  were  used:  90%  silt,  10%  kaolin;  70%  silt, 
30%  kaolin;  and  50%  silt,  50%  kaolin.   Grain  size  distribution  curves 
for  the  three  soil  mixtures  are  plotted  on  Figure  1.  Listed  in  Table  5 
are  the  Atterberg  limits  and  AASHTO  classification  for  the  soil  mixes. 
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APPARATUR  AND  EXPERIMENTAL  PROCEDURES 

Soil  Mixing  and  Curing 

Tlie   selected  mix  of  clay  and  silt  was  both  dry  mixed  and  mixed  to 
the  desired  moisture  content  using  a  Patterson  Kelly  twin  shell  liquid 
solid  blender.  The  blender,  shown  in  Figure  2,  consists  of  a  rotating  V 
shell  which  keeps  the  soil  in  continual  motion.   Inside  the  shell  is  a 
high  speed  rotating  dispersion  bar  from  which  dionized  water  is  sprayed. 
Connected  to  the  dispersion  bar  is  a  graduated  cylinder  and  stopcock 
which  controls  the  amount  of  water  to  be  sprayed. 

To  use  the  blender,  h   kg  of  the  desired  air  dried  silt-clay  combin- 
ation was  first  dry  mixed  for  about  10  minutes.   Water  was  then  added  at 
a  rate  of  about  50  to  100  ml  per  minute,  after  which  the  blender  was 
stopped  and  the  sides  of  the  shell  scraped  to  remove  sticking  soil. 
Following  scraping,  the  mixer  was  run  for  7  minutes,  stopped,  and  again 
the  sides  were  scraped.  The  mixer  was  rian  for  7  additional  minutes,  after 
which  the  soil  was  sealed  inside  a  polyethylene  bag  and  cured  in  a  high 
humidity  barrel  for  two  days.  When  adding  water  to  the  mix  at  hygroscopic 
water  content,  experimentation  indicated  that  adding  0.^%   less  v;ater  than 
calculated  gave  within  0.25^  of  the  desired  moisture  content. 

It  should  be  noted  that  enough  soil  was  mixed  at  one  time  so  that 
two  samples  could  be  compacted.  Thus  one  sample  co^ild  be  used  for  the 
frost  heave  test  and  the  other  cut  up  for  freeze  drying. 
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FIGl.«E  2   PATTERSON  KELLY  BLENDER 
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Compaction 

Description  of  Compactor 

As   shown   in  Figure   3,   an  electrically   driven  semi-automatic  knead- 
ing compactor  made  by   the  August  Manufacturing  Company  was   used.     As 
mentioned  by  Bailey   (1976),   this   type   of  compactor  seems   to  simulate 
actual  field  compaction   conditions.      Also,   since  it   is   mechanized,    the 
compaction   foot   gives   very   repetitious    foot   pressures,   and  hopefully   a 
very  uniform  sample.      This   is  highly   desirable  where  very  small  sajnples 
are  to  be  taken   from  the   compacted  soil   for  pore  size   distribution 
studies . 

Coil  samples  were   compacted  by   a  triangular  shaped   foot  with   the 
sample   rotated  60°  between  application.     The   foot  pressure  was   varied 
by   a  pneumatic-hydraulic   control   system  using  a  standard  air   regulator 
and  pressure   gauge.      At   the  present   time,    the   actual   foot  pressure   versus 
gauge  pressure   is   somewhat   difficult  to  accurately  quantify.      However, 
the  maximum  dynamic   foot  pressure  versus   gavige  pressure   could  be   deter- 
mined using  a  proving  ring.     This   relation  is   plotted  in  Figure   h.      oince 
this   calibration   varies  with   time   and  use  of  the   compactor,   the   reported 
foot  pressures   should  be   viewed  as   approximate. 

Compaction  Procedure 
Tliree  nominal  energy  levels   of  compaction  were  used  to  give   a  range 
of  soil   densitjr   conditions.      The   gauge   readings    for  the   compactive  levels 
are   as   follows : 

1)  3.5   and  i+0  psi    for  90;^   silt  -  10^  kaolin 

2)  h,   8.5   and  ko  psi   for  10%  silt   -   30f»  kaolin 

3)  h   and  8.5  psi   for  50%  silt  -   '?0%  kaolin. 
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FIGURE    3     AUGUST  MANUFACTURING  COMPANY 
KNEADING  COMPACTOR 
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The  middle  range  (8.5  psi)  of  compactlve  effort  was  selected  so  that 
the  maximum  dry  density  for  the  90%  silt  -  10%  kaolin  mixture  closely 
approximated  that  obtained  by  an  impact  hammer  using  AASHTO  test  pro- 
cedure T99-70,  method  "A",  except  that  fresh  soil  samples  were  used  for 
each  point. 

As  mentioned  by  Bailey  (1976),  there  is  no  standard  method  for  this 
type  of  compaction  with  soil.  Therefore,  the  following  procedure  was 
used.  The  sample  was  compacted  in  a  Lucite  mold  4.0  in.  in  diameter  by 
4.58  in.  high,  which  was  supported  in  a  ring  as  shown  in  Figure  5.  The 
soil  sample  was  compacted  in  5  equal  layers  with  30  blows  per  layer. 
Before  compaction,  each  layer  was  lightly  tamped  with  a  piston.  After 
the  layer  was  compacted,  the  top  of  it  was  scarified  before  compacting 
the  next  layer.   If  this  was  not  done,  cracking  tended  to  appear  between 
layers  for  soil  compacted  on  the  wet  side  of  optimum.   Once  compaction 
was  completed,  the  sample  and  mold  were  removed  from  the  compaction 
holder  and  the  soil  was  ready  for  saturation. 

Saturation 

To  maximize  the  frost  heave,  the  freezing  samples  were  brought  as 
close  to  complete  saturation  as  practicable  while  still  retaining  a 
constant  volume  condition.  To  accomplish  this,  the  saturation  device 
shown  in  Figures  6  and  7  was  used. 

To  saturate  the  sample  enclosed  in  the  Lucite  tube,  caps  containing 
porous  stones  and  qulcklock  connections  were  placed  on  each  end.   The 
assembly  was  then  connected  for  two  days  to  an  inflow  and  outflow 
reservoir  of  the  saturation  device.  The  device  was  designed  so  that 
both  reservoirs  could  be  independently  pressurized  and  also  so  that  a 
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FIGURE    5     COMPACTION  MOLD.  HOLDER.  AND  RSTON 
TAfi^PER 
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FIGURE  6    COJ^PACTION  MOLD  AND  END  CAPS  USED  FOR 
SATURATION 
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FIGURE   7      SATURATION  DEVICE 
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vacutun  could  be  connected  to  the  outflow  reservoir. 

The  usual  procedure  was  to  apply  a  pressure  to  the  inflow  reservoir 
so  that  a  gradient  would  exist  in  the  sample.  In  some  cases,  the  inflow 
and  outflow  reservoirs  were  both  pressurized,  so  that  air  bubbles  in  the 
sample  would  be  compressed  and  hopefully  flow  out  of  the  sample.  The 
gradient  used  for  saturation  depended  on  the  soil,  but  in  no  case  exceeded 
about  5  psi.  For  some  soils,  it  was  found  that  applying  a  vacuum  at  the 
top  of  the  sample  also  aided  in  saturation. 

For  almost  all  the  soil  mixes  it  was  possible  to  obtain  saturation 
of  95%  or  better.   Some  problems  were  encountered  with  the  90%  silt  - 
10%  kaolin  mixture  compacted  on  the  wet  side.  However,  soil  mixtures  of 
70%  silt  -  30%  kaolin  and  50%  silt  -  50%  kaolin  compacted  at  optimum  and 
the  wet  side  had  initial  compacted  saturations  of  close  to  95%,  and  thus 
no  problems  were  encountered.   Soil  mixes  compacted  on  the  dry  side 
seemed  to  saturate  well  by  using  a  vacuum. 

Frost  Heave  Test 

At  the  present  time  there  is  no  universal,  simple,  frost  heave  test. 
The  U.  S.  Army  Cold  Regions  Research  and  Engineering  Laboratory  (CRREL) 
test  seems  to  be  the  most  universal,  however  it  is  very  time  consuming. 
The  test  is  based  on  a  slowly  descending  frost  line,  controlled  by 
constantly  adjusting  the  room  temperature. 

A  more  practical  test  is  the  rapid  freeze  test  developed  by  Kaplar 
(1968)  and  Zoller  (1973) .  This  test  uses  a  constant  room  temperature  to 
freeze  the  sample  and  only  takes  2  days.  A  problem  with  th6  setup  that 
these  gentlemen  used  is  that  the  sample  was  compacted  and  frozen  in 
Lucite  rings.  The  rings  were  selected  to  reduce  side  wall  friction. 
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however,  they  cause  problems  in  saturating  soils  of  low  permeability. 
Since  some  of  the  soil  mixes  used  in  this  study  are  of  low  permeability, 
different  sample  containers  were  designed. 

The  sample  container,  shown  in  Figures  8  and  9,  consists  of  a  17  in. 
by  17  in.  by  6.75  in.  block  of  Dow  Styrofoam  in  which  two  4  in.  diameter 
cylinders  were  bored  to  hold  the  samples.   In  order  to  prevent  friction 
between  the  cylinder  and  S£imple,  the  sides  were  lined  with  teflon  tape. 
Attached  to  the  top  of  the  Styrofoam  were  two  brackets  which  held  dial 
gauges  to  measure  sample  heave.  Also,  porous  stones  were  mounted  in  the 
cylinder,  so  that  the  top  of  the  soil  sample  protruded  slightly  above 
the  top  of  the  Styrofoam  block.   Styrofoam  was  used  as  the  container 
because  of  its  low  thermal  conductivity,  and  thus  it  was  felt  that  the 
sample  would  be  less  likely  to  freeze  to  the  side  of  the  cylinder,  as  is 
possibly  the  case  when  using  a  Lucite  cylinder. 

The  freezing  container  with  samples  enclosed,  was  placed  inside  a 
Styrofoam  lined  box  which  contained  deionized  water  at  4  C.  The  height 
of  the  sample  container  was  designed  so  that  it  would  float  in  the  water 
bath  with  a  water  level  of  approximately  0.75  in.  above  the  bottom  of 
the  sample.  Floating  the  container  gave  a  very  simple  and  constant 
regulated  water  supply.  To  maintain  the  water  temperature  of  4  C,  a 
temperature  control  bath  and  circulator  were  used.   The  complete  setup 
is  shown  in  Figure  10. 

The  entire  assembly  was  placed  inside  a  cold  room  at  -6  C  for 
freezing.   The  average  temperature  of  the  cold  room  was  monitored  by 
placing  thermocouples  inside  a  glycerol  filled  bottle  which  rested  on 
the  sample  container.   The  water  bath  temperature  was  also  monitored 
using  thermocouples.   To  read  the  thermocouples  a  digital  voltmeter  and 


-i  ''■.:> 


39 


FIGURE  8     FROST  HEAVE  SAMPLE  CONTAINER  WITH  DIAL 

GAUGES  MOUNTED 
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FIGURE    9     POROUS  STOKES  MOUMTED  IN  FROST 
HEAVE  CONTAINER 
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FIGURE     10      CONTAirsJER   AND  SAA5PLES  DURING 
FROST  HEAVE  TEST 
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Oraega  Engineering  Inc.  cold  junction  coropensator  were  used.   The 
compensator  serves  as  an  electronic  zero  degree  reference,  thus  taking 
the  place  of  an  ice  bath. 

Testing  Procedure 

After  completion  of  saturation,  the  sample  and  end  caps  were  dis- 
connected from  the  saturation  device  and  placed  inside  a  refrigerator  set 
at  4  C.  After  tempering  overnight,  the  end  caps  were  taken  off  the  sample 
and  the  soil  was  extruded  from  the  Lucite  mold.   The  sample  was  then 
pushed  inside  the  Styrofoam  cylinder.   The  cylinder  had  essentially  the 
same  diameter  as  the  sample,  however,  very  slight  pressure  was  needed  to 
insert  the  sample.   To  aid  in  insertion,  the  teflon  cylinder  liner  was 
lightly  coated  with  silicon  grease. 

Following  insertion  of  the  soil,  dial  gauges  were  mounted  above  the 
sample  and  freezing  was  begun.   Freezing  at  the  top  of  the  sample  was 
initiated  by  leaving  the  soil  exposed  for  one  hour.  After  an  hour  had 
passed,  a  brass  cover  plate  along  with  a  surcharge  load  of  0.5  psi  to 
simulate  pavement  weight  was  placed  on  top  of  the  sample.   The  surcharge 
weights  were  mounted  in  an  elevated  position  to  allow  better  thermal 
transfer  between  the  soil,  brass  plate,  and  ambient  temperature  of  -6  C. 
On  initiation  of  freezing,  heave  readings  were  normally  noted  at  about 
4  hour  increments. 

After  freezing  for  2  days,  the  sample  and  container  were  removed 
from  the  cold  room.   The  soil  was  removed  from  the  container,  with 
photographs  of  the  soil  and  ice  lenses  and  water  contents  then  being 
taken.   Following  this,  two  new  samples  were  placed  in  the  Styrofoam 
container  and  the  process  was  repeated. 
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Freeze  Drying 

The  Process 

For  effective  measurement  of  pore  size  distribution  using  mercury 
intrusion  the  soil  must  be  almost  completely  dry.   If  compacted  soil  is 
dried  by  air  or  oven>  surface  tension  forces  produced  by  the  air-water 
menisci  generally  causes  large  amounts  of  shrinkage  and  thus  structural 
disturbance.  To  avoid  this  problem,  critical  region  drying  or  freeze 
drying  may  be  used. 

The  critical  region  process,  as  used  by  Bhasin  (1975),  subjects  the 
soil  sample  to  increasing  pressure  and  temperature  such  that  the  moisture 
remains  in  a  liquid  state  until  it  passes  the  critical  point  of  water, 
where  it  is  transformed  into  a  vapor.   Since  no  phase  interface  occurs 
within  the  critical  region,  there  are  no  surface  tension  forces  to  cause 
disturbance  to  the  sample.  For  the  present  study,  this  method  was 
initially  used,  however  destruction  of  the  soil  structure  resulted  due 
to  a  reaction  between  dolomite  in  the  silt  and  kaolinite.   See  Appendix 
B  for  further  details. 

With  freeze  drying,  the  soil  sample  is  quickly  frozen  in  liquid 
nitrogen  and  then  placed  in  a  vacuum  where  moisture  in  the  form  of  ice 
is  removed  by  sublimation.   Since  there  are  no  surface  tension  forces 
from  air-water  menisci,  the  sample  should  display  little  shrinkage  by 
this  method  of  drying.   Freeze  drying  was  successfully  used  by  Ahmed, 
Lovell  and  Diamond  (1974)  and  has  been  simplified  by  Zimmie  and  Almaleh 
(1976).  For  this  study  the  procedure  used  by  Zimmie  et  al.  seemed  to 
work  well. 
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Apparatus 
As  shown  in  Figures  11  and  12,  the  basic  components  of  the  freeze 
drying  apparatus  consist  of: 

1)  Wire  cage  to  hold  samples  during  freezing  and  sublimation. 

2)  Large  Dewar  flask  containing  liquid  nitrogen  into  which  the 
sample  and  cage  are  dipped. 

3)  Desiccator  '.«.th  inside  wire  support  from  which  sample  cage 
is  suspended  during  sublimation. 

4)  Vacuum  pump,  capable  of  evacuating  to  less  than  0.01  mm  of 
mercury . 

5)  Condenser,  which  is  hooked  between  the  desiccator  and  vacuum 
pump  in  order  to  prevent  vapor  from  the  soil  entering  the 
pump.   Liquid  nitrogen  contained  in  a  Dewar  flask  was  used  to 
cool  the  condenser. 

Procedure 
Samples  used  for  freeze  drying  were  obtained  by  quartering  the 
compacted  soil  with  a  knife  and  then  carefully  trimming  a  cube  of  soil 
with  a  razor  blade.   It  was  felt  that  a  razor  blade  produced  less 
disturbance  than  the  small  tube  sampler  used  by  Ahmed  (1971) .  After 
trimming  the  sample,  the  cubes  were  placed  in  the  cage  and  dipped  into 
liquid  nitrogen  for  about  5  minutes.   It  was  found  that  to  freeze  dry 
the  cubes  successfully,  they  could  be  no  larger  than  about  8  mm  square. 
Larger  samples  often  cracked  during  freezing,  and  displayed  greater 

shrinkage.  The  small  size  was  not  a  handicap  since  samples  -to  be  used 

3 
in  the  mercury  porosimeter  could  be  no  larger  than  about  0.4  cm  in 

volume. 
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FIGURE   ir       SAMPLE   CAGE  SUSPENDED  IN  DESICCATOR 
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FIGURE  12     FREEZE  DRYING  APPARATUS  WITH  VACUUM 
PUMP,  CONDENSER  AND  DESICCATOR 
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After  freezing,  the  sample  and  cage  were  quickly  placed  in  the 
desiccator,  which  was  then  evacuated  by  a  vacuum  pump  and  sublimation 
had  thu3  begun.  The  vacuum  pump  was  left  running  for  10  hours  with  the 
Dewar  flask  for  the  condenser  being  periodically  refilled  with  liquid 
nitrogen.   Zimmie  et  al.  (1976)  found  for  kaolinitic  clay  that  only 
about  5  to  6  hours  were  needed  for  sublimation.   However,  to  be 
conservative,  10  hours  was  used  in  this  study. 

The  basic  difference  between  the  method  used  by  Zimmie  et  al.  and 
Ahmed  was  that  Ahmed  kept  his  sample  container  cooled  to  below  0  C 
during  sublimation.   However,  as  Zimmie  et  al.  point  out  "...  sublimation 
is  an  evaporative  process  accompanied  by  cooling,  so  although  the 
specimen  container  is  exposed  to  room  temperature,  the  sample  temperature 
will  remain  well  below  freezing".   It  is  important  though,  that  the 
sample  not  touch  the  glass  desiccator  since  heat  conducted  through  the 
glass  could  cause  melting  of  the  sample.   For  this  reason,  the  sample 
cage  was  suspended  in  the  desiccator.   The  advantage  of  not  cooling  the 
sample  during  sublimation  is  that  sublimation  will  proceed  at  a  much 
faster  rate  at  higher  temperatures. 

Once  sublimation  had  been  completed,  the  sample  was  removed,  placed 
in  a  labeled  glass  jar  and  stored  in  a  desiccator  containing  anhydrous 
magnesium  perchlorate,  so  as  to  remove  any  remaining  moisture. 

The  degree  of  shrinkage  of  freeze  dried  samples  could  be  determined 
by  mercury  displacement  in  the  porosimeter  device.   For  the  soils 
studied,  it  was  found  that  not  more  than  3%  volume  change  resulted  from 
the  freeze  drying  process. 
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Pore  Size  Determination 

Concept 
Pore  size  distribution  determination  using  mercury  intrusion  is 
based  on  the  concept  that  a  non-wetting  liquid,  such  as  mercury  with 
soil,  enters  an  empty  capillary  only  under  pressure,  with  Increasing 
pressures  needed  to  intrude  smaller  pores.   If  cylindrical  pore  shape 
is  assumed,  then  knowing  the  contact  angle  between  soil  and  mercury  and 
the  surface  tension  of  mercury,  the  absolute  pressure  of  intrusion  can 
be  used  to  calculate  pore  diameter  using  the  Washburn  (1921)  equation: 

-  4  o  cos  6 
d (9) 

where  d  =  pore  diameter  in  cm, 

a  =  surface  tension  of  mercury  in  dynes/cm, 
9  =  contact  angle  between  mercury „and  soil, 
p  =  absolute  pressure  in  dynes/cm  . 

Diamond  (1970)  measured  the  contact  angle  between  mercury  and 
kaolinitic  and  illitic  soils  by  the  sessile  drop  method  and  found  it  to 
be  within  a  degree  of  an  average  value  of  147  degrees.   For  the  surface 
tension  of  mercury,  the  value  found  by  Kembell  (1946)  of  484  dynes/cm 
at  25  C  was  used  for  this  study. 

Based  on  the  computed  pore  size  diameter  and  the  amount  of  mercury 
intruded  at  the  desired  pressure,  a  pore  size  distribution  can  then  be 
calculated. 

Apparatus 
The  apparatus  used  in  mercury  pore  size  determination  consists  of 
a  penetrometer;  filling  device,  vacuum  pump,  McLeod  gauge,  mercury 


49 


manometer  and  control  board;  and  an  Aminco  porosiraeter. 

llie  glass  penetrometer  consists  of  a  bulb  to  liold  the  sample  and  a 
graduated  stem.  The  stem  is  marked  in  increments  of  0.002  ml  f>o  as  to 
measure  a  total  of  0.200  ml  intrusion  of  mercury  into  the  sample.  The 
bulb  of  the  penetrometer  is  sealed  using  a  stainless  steel  end  cap  and 
teflon  locking  ring.  Figure  13  illustrates  the  assembled  penetrometer 
and  separate  pcirts. 

Shown  in  Figure  14  is  the  filling  device,  McLeod  gauge  and  control 
board.   Tlie  filling  device  consists  of  a  two  chambered  glass  tube  in  which 
the  penetrometer  is  inserted  and  sealed.   Connected  to  the  device  is  a 
control  board  which  is  also  connected  to  a  mercury  manometer,  McLeod 
gauge  and  vacuum  pump.   The  mercury  manometer  is  used  to  measure  vacuum 
pressures  from  atmospheric  down  to  1  mm  of  mercury,  whereas  the  McLeod 
gauge  is  used  to  measure  pressures  below  1  mm  of  mercury.   Also  fastened 
to  the  filling  device  is  a  side  arm  containing  mercury,  along  with  a 
reserve  mercury  reservoir  connected  via  a  teflon  stopcock. 

The  final  piece  of  equipment  is  an  Aminco  poro£;imeter,  which 
measures  mercury  intrusion  at  pressures  greater  than  atmospheric. 
Pressure  for  intrusion  is  generated  by  the  porosimetcr  using  an 
electrically  driven  hydraulic  pump,  and  is  measured  by  two  Bourdon 
pressure  gauges  of  1000  and  15,000  psi  capacity.   Tlie  amount  of 
intrusion  is  measured  by  a  stainless  steel  needle  which  follov/s  the 
position  of  the  mercury  meniscus  in  the  capillary  stem  of  the  penetrometer. 

Test  Procedure 
The  test  procedure  used  in  this  study  was  similar  to  that  used  by 
Ahmed  (1971)  and  Bhasin  (1975).  Therefore,  the  method  is  only  briefly 
described  here. 
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FIGURE   13      PENETROMETER   DISASSEMBLED  AND  ASSEMBLED 
WITH  SAMPLE 
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FIGURE  14    FILLING  DEVICE.  McLEOD  GUAGE  AND  CONTROL 
BOARD 
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A  sample  to  be  intruded  v/as  first  trimmed  with  a  razor  blade  so  as 
to  have  a  void  volume  less  than  0.200  ml.  Following  trimming,  the  sample 
was  weighed,  placed  in  the  penetrometer,  and  then  weighed  again.  The 
penetrometer  was  inserted  in  the  filling  device  and  the  complete  assembly 
was  evacuated  to  between  0.01  and  0.02  mm  of  mercury,  which  normally 
took  about  20  minutes.   Once  filling  pressure  had  been  reached,  the  arm 
of  the  device  was  raised  such  that  mercury  completely  covered  the 
penetrometer.   The  device  was  then  shut  off  from  the  vacuum,  and  air  was 
let  in  to  raise  the  pressure  to  20  to  40  mm  of  mercury,  which  resulted 
in  mercury  flowing  into  and  filling  the  penetrometer.  After  filling  the 
penetrometer  the  arm  was  lowered,  and  pressure  was  incrementally 
increased  in  the  device  with  the  mercury  meniscus  noted  for  each 
pressure. 

Upon  reaching  atmospheric  pressure,  the  penetrometer  was  removed 
from  the  filling  device,  weighed,  and  placed  stem  down  inside  the 
porosimeter.   The  pressure  of  the  porosimeter  was  then  incrementally 
raised  until  it  reached  its  maximum  of  between  13,900  and  14,300  psi, 
with  penetrometer  probe  readings  noted  for  each  pressure  increment. 
It  should  be  noted  that  the  weight  of  the  mercury  in  a  "stem  down" 
position  caused  the  sample  to  be  at  less  than  atmospheric  at  initiation 
of  the  test.   Therefore,  to  start  the  test  at  atmospheric,  the  pressure 
in  the  porosimeter  had  to  be  raised  to  about  4  psi.  Following  completion 
of  the  porosimeter  run,  the  mercury  was  removed  from  the  penetrometer 
and  both  were  cleaned  before  reusing. 
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RESULTS  AND  DISCUSSION  OF  RESULTS 

Compaction 

As  mentioned  previously,  three  levels  of  kneading  compactive  effort 
were  used.  The  compaction  efforts  were  selected  such  that,  for  the  90% 
silt  -  10%  kaolin  mixture,  the  maximum  dry  density  for  the  middle  effort 
(8.5  psi  gauge  pressure),  closely  approximated  that  obtained  in  AASHTO 
test  procedure  T99-70,  Method  "A".  The  preliminary  compaction  curves 
for  each  of  the  three  soil  mixtures  are  shown  in  Figures  15,  16  and  17. 
It  should  be  noted  that  the  compactive  effort  referred  to  in  the 
Figures  and  throughout  the  study  are  gauge  pressures  for  the  kneading 
compactor. 

During  the  course  of  the  study  there  was  sometimes  a  slight 
difference  between  the  preliminary  compaction  curves  and  the  water 
content-density  values  obtained  during  the  actual  frost  heave  test 
program.   This  difference  is  probably  due  to  the  varying  room 
temperature  and  humidity  conditions  at  the  time  of  compaction.  The 
changing  room  conditions  were  a  result  of  the  preliminary  curves  being 
obtained  during  the  winter,  whereas  some  of  the  freezing  test  samples 
were  compacted  in  the  early  spring.  Listed  in  Table  6  are  the  water 
contents  and  dry  densities  for  the  samples  tested  during  the  freezing 
program. 

Although  sample  densities  for  similar  water  contents  varied 
slightly  over  the  course  of  the  study,  there  was  little  variance  in 
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Table  6.  Compaction  Parameters  For  Soil  Samples  Tested 
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ample 

:  %Silt-XKaolin 

Cotnpactive 

Water 

Dry  Density 

Dry  Density 

umbez 

Effort 

Content 

After 

(psl  Gauge 

Saturation 

Pressure) 

(%) 

(pcf) 

(pcf) 

1 

90-10 

8.5 

18.0 

108.7 

110.8 

2 

90-10 

8.5 

17.1 

110.3 

- 

3 

90-10 

8.5 

14.6 

110.8 

- 

4 

90-10 

8.5 

14.6 

111.0 

- 

5 

90-10 

8.5 

11.5 

108.6 

- 

6 

90-10 

8.5 

11.2 

108.5 

— 

7 

90-10 

8.5 

10.5 

106.9 

- 

8 

90-10 

40 

14.0 

114.5 

- 

9 

90-10 

40 

12.1 

120.1 

- 

10 

90-10 

40 

8.6 

118.0 

- 

11 

90-10 

40 

8.2 

117.2 

- 

12 

70-30 

4.0 

18.7 

110.8 

113.5 

13 

70-30 

4.0 

17.5 

112.9 

114.5 

14 

70-30 

4.0 

17.6 

112.5 

113.4 

15 

70-30 

4.0 

15.3 

111.4 

- 

16 

70-30 

4.0 

14.8 

103.3 

- 

17 

70-30 

8.5 

18.1 

112.7 

- 

18 

70-30 

8.5 

17.4 

113.4 

- 

19 

70-30 

8.5 

15.7 

115.4 

117.6 

20 

70-30 

8.5 

15.8 

116.1 

- 

21 

70-30 

8.5 

13.4 

115.0 

- 

22 

70-30 

8.5 

13.7 

114.4 

- 

23 

70-30 

40 

14.1 

119.4 

- 

24 

70-30 

40 

13.1 

121.2 

- 

25 

70-30 

40 

11.9 

120.0 

- 

26 

50-50 

4.0 

21.2 

104.8 

107.1 

27 

50-50 

4.0 

19.1 

106.8 

- 

28 

50-50 

4.0 

18.3 

103.6 

- 

29 

50-50 

8.5 

20.4 

108.4 

- 

30 

50-50 

8.5 

20.3 

107.3 

108.1 

31 

50-50 

8.5 

18.2 

111.4 

- 

32 

50-50 

8.5 

18.0 

111.5 

- 

33 

50-50 

8.5 

16.7 

104.7 

- 

34 

50-50 

8.5 

17.4 

106.1 

- 
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density  between  the  freeze  drying  and  frost  heave  samples  that  were 

compacted  at  the  same  time.   In  almost  all  cases  there  was  less  than 

3 

0.5  lbs/ft  difference  between  the  two  samples,  with  typical  differences 

3 
of  0.2  lbs/ft  or  less.  The  closeness  of  the  freeze  drying  and  frost 

heave  sample  densities  is  due  to  the  samples  coming  from  soil  mixed 

and  cured  in  a  single  batch,  and  also  because  of  the  high  degree  of 

replication  of  the  kneading  compactor. 


Saturation 

As  listed  in  Table  7,  almost  all  soil  samples  were  saturated  to  a 
degree  of  saturation  of  95%  or  better.   The  only  soil  samples  which  could 
not  be  brought  to  a  high  degree  of  saturation  were  the  907,   silt  -  10% 
kaolin  soil  mixtures  at  wet  of  optimum,  using  8.5  psi  compactive  effort, 
and  wet  of  optimum  and  optimum  for  the  40  psi  compactive  effort. 

Table  6  shows  that  seven  of  the  samples  displayed  higher  densities 
at  the  end  of  saturation.  This  was  due  to  the  sample  consolidating 
slightly  when  saturation  pressures  were  applied.  Determination  of  any 
change  of  volume  after  saturation  was  accomplished  by  trimming  two 
sections  of  the  freeze  drying  sample.   These  sections  were  weighed, 
coated  with  wax,  and  then  weighed  both  in  air  and  submerged  in  water, 
so  that  the  sample  volume,  density,  and  degree  of  saturation  could  be 
calculated. 

Freeze  Drying 
Table  7  lists  the  void  ratios  for  samples  both  before  and  after 
freeze  drying.   The  post-freeze  drying  void  ratios  were  determined 
during  the  mercury  intrusion  runs,  based  on  the  amount  of  mercury 


-1:' 


:,'        ■  ;^-.;    .' 


•^       -.11         C' 


■9v!£.t 


59 


Table  7. 

Compaction 

and  Frost  Heave  Parameters 

For  Soil 

Samples  Tested 

Sample 

Saturation 

Void  Ratio 

Void  Ratio 

Volume 

24  Hour 

Number 

After 
Saturation 

After  Freeze 
Drying 

Change 

Frost  Heave 

(Z) 

«) 

(mm) 

1 

89 

0.53 

0.54 

0.8 

12.2 

2 

95 

0.54 

0.51 

-2.0 

11.2 

3 

96 

0.53 

0.54 

0.2 

31.8 

4 

99 

0.53 

0.56 

2.2 

33.5 

5 

99 

0.56 

0.59 

1.7 

39.4 

6 

95 

0.57 

0.60 

2.0 

30.0 

7 

99 

0.59 

0.63 

2.3 

46.7 

8 

84 

0.48 

0.43 

-3.8 

8.4 

9 

84 

0.41 

0.44 

1.8 

11.3 

10 

99 

0.44 

0.45 

0.7 

33.0 

11 

95 

0.45 

0.47 

1.5 

24.6 

12 

95 

0.49 

0.47 

-1.3 

10.7 

13 

97 

0.47 

0.42 

-3.7 

9.4 

14 

93 

0.49 

0.48 

-0.8 

9.0 

15 

98 

0.52 

0.49 

-2.1 

34.8 

16 

? 

0.63 

0.63 

-0.3 

45.7 

17 

98 

0.50 

0.46 

-2.6 

9.3 

18 

96 

0.49 

0.49 

0.3 

6.9 

19 

94 

0.44 

0.42 

-1.5 

6.6 

20 

95 

0.45 

0.43 

-1,5 

5.7 

21 

100 

0.47 

0.46 

-0.8 

23.6 

22 

97 

0.48 

0.49 

1.3 

16.8 

23 

94 

0.41 

0.43 

0.8 

5.1 

24 

94 

0.39 

0.38 

-1.0 

5.1 

25 

89 

0.41 

0.44 

2.4 

3.2 

26 

95 

0.57 

0.57 

0.4 

8.5 

27 

96 

0.57 

0.53 

-2,6 

6.7 

28 

95 

0.62 

0.62 

-0.1 

23.6 

29 

97 

0.55 

0.51 

-2.3 

9.9 

30 

97 

0.55 

0.55 

-0.2 

5.2 

31 

98 

0.51 

0.49 

-l.A 

5.6 

32 

95 

0.50 

0.50 

-0.2 

6.1 

33 

95 

0.60 

0.56 

-2.5 

14.7 

34 

96 

0.58 

0.61 

1.7 

15.7 
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enclosing  the  sample  at  initiation  of  the  test.  This  method  assumes 
that  at  the  mercury  intrusion  filling  pressure,  no  pores  were  filled. 
For  the  soils  used,  this  assumption  was  valid.  As  can  be  seen  from  the 
Table  there  was  less  than  3%  volume  change  for  all  but  two  of  the 
samples.  These  two  samples  showed  volume  decreases  of  3.7%  and  3.8Z, 
but  these  values  may  reflect  experimental  error  rather  than  shrinkage 
due  to  freeze  drying. 

Frost  Heave 

Freezing  Results 

During  the  freezing  test  program,  frost  heave  versus  time  curves 
were  plotted  for  all  soil  samples.  Typical  heaving  curves  are  shown 
in  Figures  18  through  21.  The  curves  are  fairly  linear  for  all  soil 
combinations,  although  some  of  them  depart  slightly  from  linear  due  to 
decreasing  heave.   This  is  probably  caused  by  changing  thermal 
conditions  as  the  freezing  front  penetrated  the  sample. 

For  correlation  of  frost  heave  with  pore  size  distribution,  an 
arbitrary  heaving  value  and  criterion  were  selected.   Since  the  frost 
heave  versus  time  curves  were  particularly  linear  during  the  first  24 
hours,  the  frost  heave  after  24  hours  was  selected  as  the  evaluating  factor. 
These  values  are  shown  in  Table  7.   To  group  the  24  hour  heave  values 
into  frost  susceptible  categories,  the  classification  listed  in  Table  8, 
relative  to  the  soils  and  procedures  used,  was  chosen. 

The  effects  of  water  content  and  compaction  on  frost  heaving  rate 
are  shown  in  Figures  18  through  21,  where  the  high,  medium,  and  low 
water  contents  are  wet  of  optimum,  at  optimum,  and  dry  of  optimum. 
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Table  8.  Frost  Susceptibility  Classification  Used  for  Purposes 
of  Rating  the  Tested  Soils 


24  Hour  Frost  Heave 

Frost  Susceptibility 

(mm) 

Rating 

<  6 

Very  low 

6-11 

Low 

11  -  20 

Medium 

20  -  30 

High 

>  30 

Very  High 
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respectively.  Table  9  summarizes  the  frost  susceptibility  classifi- 
cations of  each  soil  by  the  Casagrande  criteria,  U.  S.  Army  Corps  of 
Engineers  "F"  rating,  and  the  authorb  system. 

For  the  90%  silt  -  10%  kaolin  soil  compacted  at  8.5  psl,  both  the 
dry  side  of  optimum  and  optimum  had  "very  high"  heave,  with  heave 
Increasing  somewhat  as  water  contents  decreased.  When  water  contents 
were  Increased  to  the  wet  side  of  optimum,  the  heave  decreased  to  the 
lower  "medium"  range. 

By  Increasing  the  compactive  effort  to  AO  psl,  the  heave  for  90% 
silt  -  10%  kaolin  compacted  at  optimum,  decreased  to  the  "medium"  range. 
However,  the  heave  for  the  wet  side  samples  decreased  only  slightly  and 
the  heave  for  the  dry  side  samples  were  still  in  the  "very  high"  range. 

For  the  70%  silt  -  30%  kaolin  samples,  the  dry  side  samples  had 
"very  high"  heave  for  4  psl,  decreasing  to  "medium"  to  "high"  heave  at 
8.5  psl,  and  "very  low"  at  40  psl  compactive  effort.   The  optimum  and  wet 
sample  heave  also  decreased  slightly  with  increasing  compactive  effort. 
However,  unlike  the  90%  silt  -  10%  kaolin  mixture,  the  wet  side  heave 
was  slightly  higher  than  that  at  optimum,  with  both  much  less  than  the 
dry  side  for  4  and  8.5  psl  compactive  effort. 

The  50%  silt  -  50%  kaolin  mix  behaved  much  like  the  70%  silt  -  30% 
kaolin  mix,  with  samples  compacted  on  the  dry  side  having  much  greater 
heave  than  those  compacted  at  optimum  or  wet  of  optimum.  Again,  the 
wet  side  sample  heave  was  greater  than  that  of  optimum.  Also,  Increasing 
the  compactive  effort  from  4  psl  to  8.5  psl  decreased  heave  on  the  dry 
side,  but  had  little  effect  on  samples  compacted  at  optimum  or  wet  of 
optimum. 


Table  9.   Frost  Susceptibility  Ratings  of  Soils  Tested 
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Sample   ^Sllt-^Kaolln 
Number* 


1 
2 
3 
k 

5 
6 

7 
8 

9 
10 
11 
12 
13 
lU 

15 
16 

17 
18 

19 
20 
21 
22 
23 
2U 

25 
26 
27 
28 
29 
30 
31 
32 
33 
3k 


90  -  10 
90  -  10 
90  -  10 
90  -  10 
90  -  10 
90  -  10 
90  -  10 
90  -  10 
90  -  10 
90  -  10 
90  -  10 
70  -  30 
70  -  30 
70  -  30 
70  -  30 
70  -  30 
70  -  30 
70  -  30 
70  -  30 
70  -  30 
70  -  30 
70  -  30 
70  -  30 
70  -  30 
70  -  30 
50  -  50 
50  -  50 
50  -  50 
50  -  50 
50  -  50 
50  -  50 
50  -  50 
50  -  50 
50  -  50 


Casagrande's 

Corps  of 

Author ' s 

Soil 

Engineers  Soil   Frost  Heave 

Classification 

Classification     Rating*** 

FS»« 

fu 

Medium 

FS 

Fl* 

Medium 

FS 

Fl4 

Very  High 

FS 

Fk 

Very  High 

FS 

Fk 

Very  High 

FS 

Fk 

High 

FS 

Fk 

Very  High 

FS 

Fk 

Low 

FS 

Fk 

Medium 

FS 

Fk 

Very  High 

FS 

Fit 

High 

FS 

Fl; 

Mediiim 

FS 

Fk 

Low 

FS 

Fk 

Low 

FS 

Fk 

Very  High 

FS 

Fk 

Very  High 

FS 

Fk 

Low 

FS 

Fk 

Low 

FS 

Fk 

Low 

FS 

Fk 

Very  High 

FS 

Fk 

High 

FS 

Fk 

Medium 

FS 

Fk 

Very  Low 

FS 

Fk 

Very  Low 

FS 

Fk 

Very  Low 

FS 

F3(c] 

Low 

FS 

F3(c^ 

Low 

FS 

F3{c 

High 

FS 

F3(c 

Low 

FS 

F3(c 

)       Very  Low 

FS 

F3(c 

)      Very  Low 

FS 

F3(c 

Low 

FS 

F3(c 

Medium 

FS 

F3(c 

1       Medium 

*  See  Table  6  for  compaction  variables.   Sample  numbers  correspond  to 

those  in  Tables  6  and  7. 
**  FS:  Frost  susceptible  soil 
***  The  reader  is  reminded  that  this  rating  has  not  been  applied  to 

field  frost  susceptibility.   It  applies  only  to  the  laboratory  tests 

described  herein. 
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In  sumnvary,  the  dry  side  samples  tended  to  have  greater  heave  and 
be  more  frost  susceptible.   Increasing  the  compactive  effort  did  decrease 
heave  for  soils  on  the  dry  side,  but  not  always  for  optimum  or  wet  side 
samples.  In  comparison  of  wet  side  with  optimum,  greater  heave  may  be 
displayed  for  either,  depending  on  the  soil. 

Description  of  Frozen  Soil 
As  can  be  seen  In  Figures  22  and  23,  ice  lenses  tended  to  be  much 
more  prevalent,  although  slightly  smaller,  for  the  dry  side  samples 
rather  than  optimum  and  wet  side  samples.  With  increasing  clay  contents, 
the  samples  tended  to  have  somewhat  thicker  lenses,  but  they  were  spaced 
much  further  apart.  Also,  ice  lenses  for  samples  at  optimum  or  wet  of 
optimum  and  with  high  clay  contents,  tended  to  be  more  irregularly 
oriented. 

Pore  Size  Distribution 

Determination  of  Pore  Size 

For  all  mercury  Intrusion  runs,  the  pressure  and  intrusion  data 
were  punched  on  computer  cards  and  programmed  on  the  Purdue  CDC-6500 
Dual  Mace  computer,  using  the  programs  given  in  Appendix  A.  The 
computer  tabulated  pore  size  diameter  and  cumulative  intrusion,  and  also 
plotted  limiting  pore  size  diameter,  both  with  respect  to  cumulative 
voids  per  weight,  and  cumulative  voids  per  volume. 

Drake  and  Ritter  (1945)  introduced  the  concept  of  plotting 
cumulative  voids  per  gram  versus  pore  diameter,  and  this  method  was  also 
used  by  Bhasln  (1975)  and  Ahmed  (1971).  However,  for  this  study  it  was 
felt  that  the  amount  of  pores  per  given  volume  would  be  more  relevant, 
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Dry  of  Optimum 


Optimum 


Wet  of  Optimum 


FIGURE   22       FROST  HEAVE  SAMPLES  OF   70%  SILT- 
30%  KAOLIN  COMPACTED  AT  85  PSI 
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Dry  of  Optimum 


Optimum 


Wet  of  Optimum 


FIGURE  23       FROST  HEAVE  SAMPLES  OF  50%  SILT- 
50%  KAOLIN  COf^ACTED  AT  8.5  PSI 


71 


and  the  final  analysis  was  based  on  plots  of  cumulative  voids  per  volume, 
or  "cumulative  porosity"  versus  limiting  pore  size  diameter. 

For  each  compacted  sample,  there  were  two  or  more  porosimeter  runs. 
Prior  to  the  freeze  drying  process,  two  samples  were  taken  from  the  top 
and  two  from  the  bottom  of  each  compacted  soil  sample.   It  was  found  that 
the  pore  size  distribution  curves  were  essentially  the  same,  regardless 
of  the  location  of  the  sample.  At  times  total  porosities  were  different 
in  replicate  mercury  intrusion  runs,  in  which  case  additional  runs  were 
made.  From  the  multiple  runs  for  each  compacted  sample,  a  representative 
pore  size  distribution  curve  was  selected  for  analysis. 

Comparison  of  Pore  Size  Curves 

The  pore  size  distribution  curves  are  given  in  Figures  24  through 
34,  where  the  high,  medium  and  low  water  contents  represent  wet  of 
optimum,  at  optimum,  and  dry  of  optimum,  respectively.   Comparison  of 
these  curves  shows  that  the  distribution  of  fine  pores  below  0.1  ym  was 
consistent  for  each  soil.   Significant  differences  in  distributions  of 
porosity  due  to  compactive  effort  and  water  content  did  not  begin  to 
appear  until  about  the  0.4  jam  size.   This  is  in  agreement  with 
Sridharan,  Altschaeffl  and  Diamond  (1971),  Bhasin  (1975),  Diamond 
(1971) ,  and  Ahmed  (1971) ,  who  worked  with  soils  ranging  from  Edgar 
plastic  kaolin  and  Grundite  to  Boston  blue  clay  and  Crosby  silty  clay. 

Bhasin  (1975)  suggests  that  the  fine  pores  (<  0.1  ym)  occur  either 
within  individual  particles  or  within  domains  of  roughly  face  to  face 
oriented  clay  particles.  Furthermore,  Bhasin  points  out  that  this  size 
range  conforms  roughly  to  the  intradomainal  pore  size  range  of  300A  to 
800X  as  determined  by  Diamond  (1971)  for  Edgar  plastic  kaolin.  Since 
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the  pore  size  below  0.1  \m   is  almost  in  direct  proportion  with  the  kaolin 
content  of  the  three  mixes,  one  might  conclude  that  this  pore  size  range 
is  due  to  the  clay  particles,  with  the  silt  contributing  little. 

Pore  Size  Distribution  on  Dry  Side  of  Optimum 
Figures  24  through  34  show  that  for  the  soil  samples  compacted  on 
the  di7  side  of  optimum,  there  was  a  high  percentage  of  large  pores 
(>  2  vim) ,  relative  to  samples  compacted  at  optimum  or  the  wet  side  of 
optimum.  Again,  this  is  in  agreement  with  Bhasin  (1975),  Ahmed  (1971) 
and  Diamond  (1971). 

It  is  interesting  to  note  that  for  the  90%  silt  -  10%  kaolin  soil 
mixture  when  compacted  on  the  dry  side  of  optimum  (Figures  24,  25,  26), 
a  high  percentage  of  large  pores  was  in  a  relatively  narrow  band  of 
pore  size.  For  instance,  for  the  soil  compacted  at  8.5  psi,  about  50% 
of  the  pores  were  in  a  band  between  7  and  9  ym.   This  might  suggest  a 
very  orderly  packing  between  silt  particles  and  silt-clay  aggregations 
for  the  largely  silt  soil. 

For  the  soil  mixtures  compacted  on  the  dry  side,  increasing  the 
compactive  effort  tended  to  decrease  the  amount  of  large  pores.  Again, 
the  90%  silt,  10%  kaolin  soil  showed  a  large  number  of  large  pores  in 
a  relatively  narrow  band,  however,  this  band  had  shifted  to  a  smaller 
size  in  comparison  to  the  lower  compactive  effort. 

Fore  Size  Distribution  of  Optimum  and  Wet  Side  Samples 
At  a  compactive  effort  of  8.5  psi,  the  90%  silt  -  10%. kaolin  samples 
displayed  larger  amounts  of  pore  space  between  1  to  8  ym  for  the  soil 
coi]q>acted  at  optimum  in  comparison  to  wet  of  optimum.  By  increasing  th« 
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compactive  effort  to  40  psi,  this  range  was  reduced  to  1  to  6  vun. 

In  contrast  to  the  above,  the  70%  silt  -  30%  kaolin  mixtures 
(Figures  27  through  31)  displayed  larger  amounts  of  soil  pores  in  the 
0.4  to  3  um  range  for  the  8.5  psi  compactive  effort,  and  in  the  0.4  to 
1  ym  range  for  soil  compacted  at  40  psi.   Outside  of  the  above 
mentioned  ranges,  the  amount  of  pores  was  essentially  the  same, 
whether  the  soil  was  compacted  at  optimum  or  wet  of  optimuai. 

For  the  50%  silt  -  50%  kaolin  mixture,  the  wet  side  of  optimum 
displayed  larger  amounts  of  pores  than  the  optimum  for  the  range  of 
0.1  to  3  ym  at  the  4  psi  compactive  effort  and  in  the  range  of  0.4  to 
3  ym  for  soil  compacted  at  8.5  psi.   See  Figures  32,  33  and  34. 

Relation  of  Frost  Heave  to  Pore  Size 

Regression  Analysis 
To  derive  a  prediction  equation  for  frost  heave  as  a  function  of 
pore  size,  the  method  of  linear  regression  was  used.   This  method  is 
based  on  a  model,  which  In  the  case  of  a  first  order  relation  Is 

Y^  -  e^  +  B^  x^  +  e^  (10) 

where  3  and  3-,  are  parameters  of  the  model,  X.  is  the  independent 
variable,  Y,  is  the  corresponding  observation,  or  dependent  variable, 
and  e.  is  the  departure  from  the  regression  line,  or  random  error. 

In  the  above  model,  3  >  3,  and  e.  are  unknown,  with  e,  changing 
for  each  observation  Y , .  Although  e.  is  impossible  to  predict,  3  and 
3-1  can  be  estimated  by  use  of  the  method  of  least  squares,  which  is 
based  on  minimizing  the  sum  of  squares  of  errors.  The  sum  of  squares 
of  errors  Is  equal  to 
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i=l      i=l 


To  use  the  least  squares  methods,  estimates  of  6  and  8, ,  denoted 

o      i 

by  b  and  b, ,  are  picked  so  that  when  substituted  for  3  and  B^ ,  the 

sum  of  squares  of  errors  is  minimized.   Once  the  estimates  b  and  b,  are 

o     1 

determined,  the  prediction  equation  is  written  as 


Y  =  b^  +  b^  X  (12) 

where  Y  is  the  predicted  value  of  Y. 

For  the  above  prediction  equation,  the  following  assumptions  have 

been  made. 

2 

1)  e  are  random  variables  with  mean  0  and  variance  a  . 

2)  c.  and  c  are  uncorrelated,  i  j'  j ,  so  that  cov  (e  ,  e  )  =  0. 

3)  For  the  F  test,  used  as  a  method  of  evaluation  of  the  model, 

e  is  a  normally  distributed  variable  with  mean  0  and 

2 
variance  o  . 

A  measure  of  how  well  the  prediction  equation  fits  the  actual  data 

2 
is  the  coefficient  of  determination,  R  ,  which  is  defined  as 


2  _  Sum  of  squares  due  to  regression 
Total  (corrected)  sum  of  squares 


-  2 

--^ — r-2  (13) 

£(Y  -  Y)^ 


The  coefficient  of  determination  is  usually  expressed  as  a  percentage, 

2 

100  R  ,  and  is  the  percent  variation  explained  by  the  regression  line. 
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2 

If  R  Is  0%,  then  no  variation  is  explained  by  the  prediction  equation. 

2 
On  the  other  hand,  if  R  is  100%,  then  all  variation  is  explained  by  the 

equation. 

A  second  test  for  the  applicability  of  the  model  is  the  F  test, 

which  tests  whether  3,  can  be  considered  non  zero  or  not.  For  this  test, 

the  ratio 


„  _    Mean  Square  due  to  regression 

Mean  Square  due  to  residual  variation 


is  compared  with  the  100  (1  -  a)%  point  of  the  F  {1,  (n  -  2)}  distribu- 
tion in  order  to  determine  if  B,  can  be  considered  non  zero  for  the 
available  data.   If  it  cannot,  then  the  regression  equation  has  little 
value.  A  complete  explanation  of  the  F  test  is  given  in  Draper  and 
Smith  (1966). 

For  the  actual  regression  analysis  for  this  study,  prediction 
equations  were  developed  using  the  Statistical  Package  for  the  Social 
Sciences  (SPSS)  computer  program  REGRESSION,  developed  by  Nie  et  al. 

(1975) .   The  equations  were  then  evaluated  for  their  accuracy  by 

2 

(1)  The  coefficient  of  determination  R  . 

(2)  The  overall  F  test,  with  the  test  being  significant  at 
a  =  0.5. 

(3)  The  partial  F  test  for  each  6. ,  with  the  test  being 
significant  at  a  =  0.05. 

Prediction  of  Frost  Heave 
The  pore  size  distribution  curves  shown  in  Figures  24  through  34 
were  evaluated  in  the  following  ways: 


(' 
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(1)  Calculation  of  a  mean  and  standard  deviation  for 
each  curve. 

(2)  Division  of  the  curves  into  percentage-of-intrusion 
bands.  These  bands  indicated  the  pore  size  for 
which  10%,  20%,  ...  90%  of  the  pores  were  larger. 

(3)  Division  of  the  curves  into  bands  of  pertinent  pore  size 
diameters.  From  these  bands,  the  cumulative  porosities 
were  obtained  for  pore  size  diameters  of  0.05,  0.08, 
0.1,  0.2,  0.4,  1.0,  1.5,  2.0,  3.0,  A.O  and  8.0  ym. 

Table  10  presents  the  data  for  methods  (2)  and  (3).   The  maximum  pore 
size  that  was  measured  was  300  \m. 

For  the  frost  heave  prediction  equations,  parameters  from  one  of 
the  above  three  cases  were  picked  as  the  independent  variables,  with 
frost  heave  as  the  dependent  variable.   Regression  analysis  was  then 
used  to  formulate  and  evaluate  the  prediction  equations. 

The  first  parameters  to  be  evaluated  were  the  mean  and  standard 
deviation,  since  these  values  could  be  calculated  using  the  same 
computer  program  for  tabulating  and  plotting  the  pore  size  curves. 
However,  the  two  parameters  turned  out  to  be  essentially  useless  for 
predicting  frost  heave.   It  was  noted  that  the  mean  and  standard 
deviation  could  be  greatly  altered  by  a  few  very  large  pores.   Thus, 
this  method  of  prediction  was  abandoned. 

The  second  method  of  prediction  was  to  consider  the  pore  diameters 
for  selected  percentages  of  intrusion.   This  method  resulted  in  fairly 
good  results.  For  a  one  order  linear  model,  the  best  prediction  was 
obtained  by  the  equation 
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Table  10.  Pore  Size  Distribution  Values  Used  in  Prediction  Equations 

Sample    Total     Cumulative  Porosity  for  Pore  Dinmeters  >  Indicated 

Number*  Cumulative  Diameter  (in  ym) ,  but  <  300  pm 
Porosity 

0.05  0.08  0.10  0.20  0.40  1.0  1.5  2.0  3.0  4.0  8.0 

.331  .319  .310  .285  .267  .239  .210  .179  .125  .081  .027 

.312  .295  .286  .261  .242  .203  .170  .140  .097  ,060  .027 

.322  .307  .298  .276  .259  .225  .204  .190  .169  .153  .100 

.340  .324  .313  .292  .277  .249  .235  .225  .209  .198  .110 

.352  .337  .328  .308  .290  .261  .247  .238  .224  .212  .135 

.338  .322  .313  .292  .277  .248  .226  .211  .190  .173  .110 

.359  .343  .334  .312  .294  .270  .258  .249  .236  .225  .115 

.261  .253  .244  .222  .204  .160  .121  .085  .050  .037  .025 

.280  .263  .255  .234  .213  .173  .152  .132  .110  .086  .018 

.290  .275  .266  .245  .266  .203  .196  .190  .177  .152  .011 

.294  .277  .270  .249  .235  .212  .203  .198  .182  .147  .014 

.277  .240  .217  .169  .137  .079  .046  .035  .027  .022  .018 

.251  .214  .192  .146  .121  .070  .039  .030  .027  .025  .020 

.262  .237  .216  .172  .140  .080  .049  .037  .028  .024  .018 

.267  .229  .208  .165  .137  .095  .071  .062  .055  .050  .038 

.389  .358  .338  .294  .268  .235  .209  .182  .117  .055  .020 

.270  .231  .210  .164  .139  .096  .066  .045  .026  .020  .012 

.282  .243  .222  .178  .152  .100  .069  .048  .034  .026  .017 

.243  .303  .181  .140  .112  .064  .044  .035  .027  .023  .017 

.248  .209  .189  .146  .115  .063  .046  .037  .029  .026  .017 

.270  .231  .210  .170  .141  .109  .096  .089  .080  .073  .018 

.282  .242  .222  .182  .156  .122  .107  .098  .086  .075  .017 

.243  .205  .180  .140  .111  .060  .044  .035  .029  .026  .020 

.210  .170  .148  .112  .090  .057  .043  .035  .029  .025  .019 

.261  .222  .200  .160  .136  .105  .092  .082  .044  .017  .011 

.299  .237  .205  .148  .114  .074  .060  .050  .038  .031  .022 

.278  ,218  .184  .120  .079  .051  .045  .040  .037  .033  .029 

.321  .260  .230  .180  .149  .110  .092  .078  .064  .055  .040 

.277  .212  .181  .131  .100  .070  .058  .049  .040  .032  .019 

.281  .215  .187  .136  .107  .072  .058  .050  .039  .029  .019 

.255  .180  .156  .108  .082  .057  .049  .041  .034  .030  .021 

.255  .185  .160  .110  .081  .057  .050  .044  ,039  .034  .029 

.294  .238  .209  .150  .119  .088  .075  .064  .042  .028  .008 

.315  .268  .237  .170  .130  .099  .069  .055  .037  .025  .018 


*  See  Table  6  for  Compaction  variables.  Sample  numbers  correspond  to 
those  in  Tables  6,7  and  9. 
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.340 
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.326 
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.336 
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.353 
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.365 
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.350 
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.374 
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.284 
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.303 

11 

.311 
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.301 

13 

.282 
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.298 

15 

.293 

16 

.413 

17 

.298 

18 

.309 

19 

.275 

20 

.278 

21 

.298 

22 

.310 

23 

.276 

24 

.241 

25 

.290 

26 

.339 

27 

.319 

28 

.362 

29 

.321 

30 

.325 

31 

.302 

32 

.303 

33 

.335 

34 

.353 
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Table  10. 

(Coiit' 

d) 

Sample 

Pore 

Diameter  for 

Indicated 

Number* 

Percentage  of 

Intrusion 

(in  lim) 

90% 

80% 

70% 

60% 

50% 

40% 

30% 

20% 

10% 

1 

0.113 

0.34 

1.02 

1.58 

2.17 

2.77 

3.45 

4.55 

6.90 

2 

0.087 

0.212 

0.64 

1.15 

1.63 

2.23 

2.97 

3.75 

5.60 

3 

0.087 

0.280 

0.81 

1.56 

3.00 

5.55 

7.90 

9,30 

10.5 

A 

0.091 

0.323 

1.06 

2.90 

5.70 

7.30 

8.00 

8.50 

9.00 

5 

0.101 

0.405 

1.20 

3.45 

6.30 

7.80 

8.40 

9.20 

10.3 

6 

0.097 

0.380 

1.07 

1.98 

3.80 

6.10 

8.10 

9.40 

10.8 

7 

0.096 

0.340 

1.35 

3.93 

6.10 

7.50 

8.10 

8.60 

9.60 

8 

0.073 

0.169 

0.495 

0.86 

1.21 

1.56 

1.99 

2.70 

6.80 

9 

0.082 

0.217 

0.55 

0.97 

1.63 

2.67 

4.00 

4.80 

6.10 

10 

0.083 

0.215 

0.67 

2.67 

3.95 

4.60 

4.90 

5.20 

5.60 

11 

0.07A 

0.217 

0.84 

2.82 

3.80 

4.25 

4.70 

5.10 

5.70 

12 

0.053 

0.077 

0.104 

0.150 

0.29 

0.54 

0.86 

1.25 

2.45 

13 

0.046 

0.068 

0.092 

0.126 

0.222 

0.48 

0.81 

1.17 

2.13 

14 

0.047 

0.078 

0.106 

0.169 

0.333 

0.59 

0.90 

1.28 

2.58 

15 

0.052 

0.074 

0.102 

0.151 

0.310 

0.66 

1.12 

2.27 

11.3 

16 

0.065 

0.107 

0.224 

0.73 

1.54 

2.25 

2.90 

3.43 

4.70 

17 

0.050 

0.075 

0.099 

0.142 

0.29 

0.64 

1.10 

1.62 

2.68 

18 

0.052 

0.077 

0.103 

0.168 

0.385 

0.71 

1.10 

1.63 

3.30 

19 

0.045 

0.065 

0.078 

0.117 

0.214 

0.42 

0.74 

1.18 

3.13 

20 

0.047 

0.068 

0.092 

0.127 

0.230 

0.43 

0.72 

1.12 

3.50 

21 

0.050 

0.071 

0.100 

0.160 

0.330 

0.72 

1.90 

5.60 

7.30 

22 

0.050 

0.073 

0.103 

0.180 

0.405 

0.90 

2.32 

5.00 

7.00 

23 

0.046 

0.064 

0.088 

0.120 

0.212 

0.41 

0.68 

1.10 

3.25 

2A 

0.042 

0.059 

0.077 

0.101 

0.152 

0.318 

0.66 

1.22 

4.25 

25 

0.048 

0.070 

0.096 

0.147 

0.310 

0.71 

1.72 

2.62 

3.48 

26 

0.048 

0.062 

0.077 

0.100 

0.138 

0.250 

0.53 

1.17 

3.45 

27 

0.044 

0.058 

0.074 

0.094 

0.120 

0.178 

0.290 

0.58 

5.70 

28 

0.045 

0.063 

0.083 

0.110 

0.198 

0.44 

1.02 

2.40 

10.1 

29 

0.044 

0.056 

0.074 

0.092 

0.120 

0.21 

0.45 

1.22 

4.20 

30 

0.042 

0.057 

0.071 

0.091 

0.123 

0.24 

0.50 

1.28 

3.55 

31 

0.040 

0.055 

0.066 

0.080 

0.103 

0.152 

0.31 

0.88 

3.90 

32 

0.041 

0.055 

0.066 

0.081 

0.103 

0.158 

0.31 

0.85 

6.60 

33 

0.041 

0.062 

0.081 

0.103 

0.148 

0.270 

0.68 

1.88 

3.52 

34 

0.048 

0.069 

0.093 

0.120 

0.180 

0.315 

0.66 

1.47 

3.00 
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where   Y  =  predicted  frost  heave  rate  in  ran/day, 

D,-  =  pore  diameter  such  that  40%  of  the  pores  are  larger, 

Dq-,  =  pore  diameter  such  that  80%  of  the  pores  are  larger, 
ou 

For  the  above  prediction  equation,  R  =  69.0%,  and  the  residuals (c)  are 
plotted  in  Figure  35. 

The  last,  and  most  fruitful,  method  of  analysis  was  to  consider 
the  cumulative  porosity  for  pertinent  pore  size  bands.   The  approach 
was  to  first  run  a  one  order  regression  analysis  for  each  soil,  looking 
at  all  eleven  selected  pore  size  diameters.   Based  on  these  prediction 
equations,  listed  in  Table  11,  this  approach  seemed  promising. 

The  next  step  was  to  run  regression  analyses  for  all  data 
combined.   When  considering  each  soil,  the  prediction  equation 
depended  only  on  the  large  pores.   This  is  partly  due  to  the  earlier 
mentioned  fact  that  the  amount  of  smaller  pores  was  essentially 
constant.   Therefore,  it  seems  reasonable  that  a  prediction  equation 
for  all  soils  would  also  have  to  take  account  of  the  small  pores. 
Based  on  this  premise,  various  two  order  models  were  tried,  and  the 
best  results  were  obtained  by  the  prediction  equations  listed  in 
Table  12.  The  residuals  for  these  equations  are  plotted  in  Figures 
36,  37  and  38. 

Proposal  and  Discussion  of  Frost  Heave  Prediction  Equation 

Based  on  the  results  discussed  in  the  preceding  paragraph,  the 
equation 


Y  =  -  5.46  -  29.46 


^3.0 


^0  "  ^0.4^ 


+  581.1  (X^  q)       (15) 
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where    ^  =  frost  heave  rate,  in  nnn/day, 

X_  =  cumulative  porosity  for  pores  >  3.0  ym  but  <  300  (jm  , 

X^  =  total  cumulative  porosity, 

X_  ,  "  cumulative  porosity  for  pores  >  0.4  ym  but  <  300  ym, 

is  recommended  for  predicting  frost  heave  from  pore  size  distribution. 

Possibly  a  more  "precise"  model  can  be  found  by  using  more  parameters 

however  much  can  be  said  for  simple  models. 

In  physical  explanation  of  the  above  equation,  it  is  believed  that 
the  large  pores,  quantified  by  the  value  X_  _,  offer  the  least  resistance 
for  water  to  move  to  the  freezing  front,  whereas  the  smaller  pores  rep- 
resent the  amount  of  available  energy  to  do  the  work  of  frost  heave. 
For  soils  of  the  same  type,  the  fine  pores  are  essentially  independent 
of  compaction  and  water  contents  and  thus  the  larger  pores  will  control 
the  frost  heave.   However,  when  evaluating  different  soils,  heave  will 
be  controlled  by  both  fine  and  large  pores. 

With  regards  to  the  predicted  heave  value,  the  greater  the  heave, 
the  greater  the  frost  susceptibility.   Although  there  is  no  specific 
correlation  with  field  performance  for  the  freezing  test  used,  an  indi- 
cation of  the  degree  of  frost  susceptibility  can  be  found  from  the 
prior  mentioned  author's  criterion,  shown  in  Table  8.   The  UNH  classi- 
fication, listed  in  Table  3,  is  also  somewhat  applicable  to  the  test 
results,  although  probably  quite  conservative. 

If  the  soils  studied  had  been  evaluated  solely  by  Casagrande's 
(1932)  criteria,  they  would  have  all  been  classified  as  frost  susceptible. 
When  the  soils  are  assigned  "F-ratings"  in  the  U.  S.  Army  Corps  of 
Engineers'  System  (Table  4),  they  are  classified  as  F-3  or  F-4. 
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However,  by  varying  clay  content,  water  content,  and  compactive  effort, 

frost  heave  ranged  from  "very  low"  to  "very  high".  This  range  of 
heave  can  be  predicted  by  using  the  proposed  frost  heave  prediction 
equation.  The  pore  size  based  equation  is  valuable  in  that  it  allows 
one  to  predict  the  degree  of  frost  heave  (and  relative  frost  suscept- 
ibility) as  it  is  affected  by  compaction  and  soil  variables. 
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CONCLUSIONS 

For  the  sllty  soils  studied,  the  following  conclusions  are  drawn 
with  respect  to  frost  heave,  pore  size  distribution  and  the  major 
compactive  variables. 

1)  Samples  compacted  dry  of  optimum  displayed  greater 
heave  than  those  compacted  at  optimum  or  wet  of 
optimum. 

2)  Increasing  the  compactive  effort  reduced  heave 
for  samples  compacted  dry  of  optimum. 

3)  It  was  not  possible  to  predict  whether  a  soil  will 
have  the  higher  frost  heave  compacted  at  optimum 
or  compacted  wet  of  optimum. 

4)  The  rapid  frost  heave  test,  as  used  in  these  studies, 
gives  reliable  results  with  relatively  Inexpensive 
test  equipment. 

5)  As  long  as  the  sample  to  be  dried  is  no  larger  than 
about  8  mm  square,  the  freeze  drying  process  produces 
little  volume  change  (+  3%)  of  the  sample. 

6)  For  a  particular  soil  type,  the  distribution  of  pores 
below  0.1  ym,  and  in  most  cases  below  0.4  pnit  is 
essentially  the  same  regardless  of  the  compaction 
water  content  or  effort. 

7)  The  degree  of  frost  heave  can  be  predicted  by  the  equation 
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-  5.46  -  29.46 


^3.0 


^0  ■  ^0.4 


+  581.1  (X^  q)  (15) 


where  Y    »  frost  heave  rate,  in  mm/day, 

X_  _  =  cumulative  porosity  for  pores  >  3.0  vim  but  <  300  \m, 

X«*  =  total  cumulative  porosity, 

X_  ,  =  cumulative  porosity  for  pores  >  0.4  ym  but  <  300  ym. 

8)  Frost  susceptibility  criteria  based  upon  textural 

or  grain  size  distribution  measures  have  a  more 

limited  prediction  capability  than  pore  size  ones, 

since  the  former  do  not  take  compaction  variables 

like  moisture  content  and  density  into  account. 

Such  compaction  variables  do  influence  the 

distribution  of  porosity. 


99  R 


RECOMMENDATIONS  FOR  FURTHER  RESEARCH 

The  evidence  of  this  research  supports  the  hypothesis  that  the 
distribution  of  porosity  may  be  a  superior  predictor  of  frost  heave  rate 
and  frost  susceptibility  of  fine  grained  soils.  This  type  of  work  should 
be  extended  in  the  following  categories. 

1)  Dehydration  Studies.   Both  critical  region  drying  and  freeze 
drying  should  be  experimentally  examined  over  a  wide  range  of 
soil  types  and  conditions  in  order  to  determine  where  each  can 
be  used  with  confidence. 

2)  Porosimeter  Apparatus.  The  versatility  of  the  porosimeter 
would  be  extended  if  the  equipment  were  modified  to  accom- 
modate larger  samples  of  more  coarse  material. 

3)  Freezing  Tests.  A  standard,  universal,  rapid  freezing  test 
should  be  established.  This  test  should  be  correlated  to 
appropriate  field  measures  of  frost  susceptibility.  When  this 
is  accomplished,  practical  predictions  of  field  frost 
susceptibility  can  be  determined  from  simple  laboratory  tests, 
including  pore  size  distribution. 

k)      Control  of  Frost  Susceptibility.  An  ultimate  short  range 

objective  of  experimental  research  of  this  type  would  be  the 
control  of  frost  susceptibility  of  a  given  soil  type  through 
control  of  the  compaction  variables.  This  will  require  extensive 
laboratory  and  field  testing,  but  the  potential  benefits  in 
improvement  of  the  state-of-the-art  Justify  the  effort. 
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Appendix  A 


C  THIS   hVDhWc.n    CALCUlrtltS    fiND    TAEULAl  ES   PURt:    Sll't:    DISTRieUTIDN 

C  PihRAHETER::.    hUGHij   UITH   Fl.DlTiriG    PDRE    D  i  hME  I ER   VtRoUS    CUnULHTIVE 

C  i.lTRlJSiDtl- CRhM.    CUMULhTIVE    INTPUriaM- VDLUME    fiHD    ItfCREMENTHL 

C  INTRUSinn    -VCILI-IML.       THE    NEntl   hHD    VftRInriCE    hPE    HLSG    CALCULPiTED 

PRDGRAtl   f--DRE    <;  IHPI.IT-  DUTPl.iT.  TrFF5=IHPUT.  TftPEe=DUTPUT) 

D I  riEM  V I  on    3 ■::  f  0  :■ ..  F'  k  5  O  '  <  D  ■  f' O >  »  CU  <■. 5 O >  «  SP  <  5 0 ::>  J  CI  ISP <  50  :>^  GRAPH <  4t.  1  )  I 
r:p:]RD'..-1Sl-.^HI:::Ta<;45l':>-Cv5CO 

REAH »    NDA I  A 

D3    4  00    .1=1  .HIiATA 
C  THE    ruriPD.E    OF    THI5    OT    LDDP    IS    TD    FE    APLE    TD   CALCULATE    RESULTS   EUR 

C  MORE    TMhM    one    SET    ITE    DATA    IM   A    STHOLE   RUN.       riDATA    IS    I  HE    hUMEEP 

c    nr  SETS  nr  dhTa  td  be  chlculhtetl  hdata  is  Att  itnrGER  value,  and 

C      iS  THE  FIRST  DATA  CARD  FUR  THE  SET  OF  DATA. 
C      DETERIUHhiTIDH  HE  PDRE  SIZE  DISTRIBUTIDH 

RiAD.DATE-Ih 
C      DATE  i-itiD  'AMPLE  NUMBER  CSN)  ARE  TYPED  CH  THE  FIRST  DATA  CARD  FDR 
C      EACH  CURVE. 

U::'irE  <6:.  14) 
14     fORMAT  a  Hi; 

PRiriT.    DATE-    SN 

READ.  US.  PL  ,  Dfl.  VP.  USP.  IJSPM,  RE.  PF.  SR 

c  ij:=ueight  of  sample.   ds=:pecific  gravity  of  sni  'ds.   dm^demsity  of 

C  MSRCUPY.    VP^-VDLUME    DF    PEMErPDMETER.    li:P=UEIGHT    OF    SAMPLE    AtlD   PEN- 

C  ETRDHETER.    USPM=U:-:iGHT    OF    SAMPLE    .    PEtlEIRDNE  [ER    AND   MERCURY  j 

C  PC.    =EVACUhIING    pressure.    PF=FILLING    pressure. in   mm    DF    MERCURY. 

C  SS    =STEN   REhDING    at    filling    pressure.       THE    ABOVE    VALUES    ARE    ALL    DM 

C  THE    SAr',E    DATA   CARD 

P;;'iriT.  US.  DS.  jiM.  VP.  USP.  USPMjPE.  PF  .  SR 

VS=IJS-DS 

IW=USPM    -    lISP 

VM=NM.--DM 

VSA=^VP-VM*':  1 .  +PE-'PF)-SR 

VVD-VSA  -VS 

VL1IDR=VVD--VS 

PJRnS= VVD-VSA 

PV=PE*VN 

VRI-PV'-PF 

PRir^T.KVDIJJME    SDL  IDS    IS*.  VS. -'.  ^VDLUME    VDIDS    IS^.VVD 

PrvItiT.*vnjD   RAT  ID    ISif^.  VDIDR.  •'.  ?PDRDSITY    ISs^-PDRDS 

REhD.N-M.SI 
C  N^lriTEGER    VALUE    FDP    THE    flUMBER   DF    LDU   PRES-:URE    PEADItlGS    MhDE>    NDT 

c  r.nurniriG  the  filling  pressure  rerding.     m-ihteger  value  fdr  the 

C      NUMBER  DF  HIi^H  PRESSURE  READItlGS  MADE.  flDT  CDUHI  It!G  THE  INITIAL 
C      PGRDSIMETEF'  READING.   SI  =  ir<TRUSIDN  READItIG  AT  INITIATIDN  PRESSURE 
C      IH  ML.   THE  THREE  VALUES  ARE  ALL  DN  DHE  DATA  CARD. 
C      LDU  PRESSURE  CALCULATIDNS 
tJRITE  <:f=..5) 

5  FORMAT  <'V.37X.?-'LDW  PRESSURE  CALCULATiaNS?«> 
URITE  C6-I.6.; 

6  FDRMAT  <-'--■>  S-^X.  ?!MM  DF  HG   PRESSURE   VDL  AIR    STEM  RDG   IHTRUS     I 
ANrPU-'GM   CriMUL  IN   INTRrVSA   CU  IN-'VS   DIAMETER   ;^) 

cu<n   =  0. 
CUSP a >  =  0. 

DD  100  !-?J.N+] 

READ.pa>..sa>  . 

C      PO>  AND  Sa;'  ARE  THE  PRESSURE  AND  INTRUSION  VALUES  FDR  EACH 

C      READING.  PtlJ  IS  GIVEN  IN  MM  DF  MERCURY  FOR  THE  LCW  PRtA?;URE 

C      INTRUSION  AND  PSI  FDR  THE  HIGH  PRESSURE.   SCI)  IS  IN  ML.   THE  TWO 

C      VALUES  ARE  TYPED  DN  ONE  DATA  CARD.  WITH  A  CARD  FOR  EACH  REATIIMQ. 

PS=PCI>«. 01934 

VA=RV--P<I) 

CVA=VAI-Vh 

csR=sa;'-SR 

SN=CSR-CVA 

IF  CSN.LT.O.)  GD  TU  10 

SG  =srH.-us 

GO  TD  11 

10  SG=0. 
S'S  =0. 

11  cua)  =  rua-i)  +  ^G 
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cuip";.  I  ;•  =  cuif'-i;  i- 1  >  »-.2p<  i  > 

pivmr  1  .f'.  1  .,|-.\'h.'(  t;-'.  SH.sn.cLK  !>•  :i;p<i).cn.sp<:i.>«Da> 
c  p*  I  ■-itnp'i  iiui  ppt  ;-i:upf    it(  Mh  Of  mpcupi".  p: -mipir  inti  ppessupl  in 

r  p:i.  vh- Villi  IMC  [ii    hip  in  ml.    :•!•     tpm  pi.Miiiriii.     ii  ihrkfcted 

L  KiTPU'.iDii  M]i-  tw  H  jtiipEMirir  in  t-iL.   ci"j-tDrK(- rr  II  Kripuiiuti  ,'GRAn 

I.  riP  ErtCH  iiiLpnifrir.   CU'. i '=ciJMiiLHnvE  ruppELirn  itupujinti-GPAn. 

c  :'"■  I v-^ini.r'FMt-jnHL  popligity*   cuspa  •^cumulhtive  pgpdsityj 

C  D'  iJ-LIhl  riHij   FDRE   Sli'E    DIAMETER 

1    r  iPMrtTcx,  1  on  0.4:) 

V;4 1  =VH 
100    SP=i:(;i::' 
C     HIGH  PRES'iURE  CALCUL--"!  IDNS 

URITE  Cb'?,' 
?      FaPMhiT  C.-.-j36Xj*HIGH  PRESSURE  CALCULftTIQNS»^> 

8      F'JRNniT  r:.-'  .3'':.*PPE:S-:URE   STEM  RDG   HG  CDR     CHG  HGC    IMTRUS     I 
BMTPU  GM   CUMUL  Ih   IMTR'VSh   CU  IH-'VS   DIAMETER   *> 
HC  i  =  1 1 .  * .  IJ  0 1 1  *VM  ■■  15000. 
DJ  ?00  I=M+£.Mf-h+l 
READiP<I>'Sci:) 

pa,'=p<i>+ii. 

HC-P  U  > * .  0  0 1 1  *VM •■- 1  "3  0  0 0 . 
CMC=HC-HCI 

csR=sa>-si 

SH=CSR-CHC 

cjsp(i>  =  cusp(;i-i>  +  spci> 

.S.iH--i:r)-' WS 

CMa^^-cua-D+SGH 

Da  '^<(-4.  ■♦■IS4.«i:aS'::£.^'=.':.6)«.  145)-'Pri> 

P.<:If<T    £.  P'.  1  ■ .  S'  I  ' .  HC.  CHC.  't1»  SGH.  CU''  I  >  .  SPC  I  ■'  •  CUSPC  I  )  «  D(I  > 
f:  pa  '-ItlTRUjIUfl    PPLSSIIPE    IN   PSI-    Sa'^STEM   PEhDIHG-    HC-MERCURY 

C  CnRrECTiaH.    CHC^CHhHOE    hi   MEPrUPV   CDPPErTIGH    in   Ml  .    SM^ilUPRECTED 

C  IHTPI.ISIDH    m    ilL.    SGH-CUPPECTED    IHTPUSIUM    -'GRAM.    CUa  ::'^lUMULATIVE 

C  rriPPb'LTEO    irMKU';II]ri.-GRAM.       SP a  )^- 1 hCPEMEMT;^L    POROSITY.    nUSP';i>  = 

C  CUMULATIVE    POROSITY) 

£  FORMhT    <;MuP1Ci.4) 

riCI=HC 

£00      si=sa) 

ri.:  300  i-^£.iHMt-i 

3 1 1 M         C  (.  I  J  -  ALDG 1  0 '.:  I'l  <  1  '  > 

C      THE  PURPOSE  OP  THIS  DO  LOOP  IS  TO  CnNVERT  THE  PORE  SIZE  DIAMETER 

C      TO  THE  COG  OF  THE  DIAMETFR 

C'CL  PLDTl  '..0.4.  10.  5.  SO) 

CALL  PLO  T  £  ■:  GRAPH .  45 1 .  3 .  .  -£ .  .  .  4 .  0 .  > 

CALL  PLOTS  a H* . C . CU . H+M > 

3  FGRMhT  aH1.35X.54l-IPL0T  OF  LIMITIMG  PORE  DIAMETER  VS  CUMULATIVE  IM 
9TRUSI0M--.-:) 

USITEtA.S' 

CALL  PC0T4  '::5.  IMTRU) 
URITE  <6.£0> 
URITE<;fc..4> 
£0     FORMAT  '::iOX.i»;0.  01^.  17X.»sO.  1?*.  18X.?!1.  LV.  1?X.'*10.;'.  1?X.  ?«1  00.  ?«>  16X 
C*i  0000.  jt'.i 

4  FORMAT  '::iH0.4'=.X.33HLIMITinG  PORE  DIAMETER  IN  MICRDriS.> 
C:'-;LL  PLOT  1  <  0.  4.  t  0.  5>  £0;' 

CALL  PLUJd    '::PaRO.  451 .  3.  .-£...  4.  0.  ) 
C-iLL  PLDT3  ■:;  1 H* .  C ,  CUSP .  tHM > 
IS     FORMAT  aHl.£':.X.7£HPLDr  OF  LIMITItIG  PORE  DIAMETER  VS  CUMULATIVE  IN 
DT^'USIOH.  VOCUHE  DF  SAMPLE/-') 
URITE  ^A. 1£) 
CALL  PLDT4  CS- INTRU) 
URITE  ^6.£0) 
WR I  TEC 6. 4) 

CALL  PLOT  1  '::0.4.  10.5»£0) 
CALL  PLDT£  •::HISTD.  451 .  3.  > -2.  j ,  2..  0.  > 
CALL  PLOTS  ':  1M*.C.SP.M-:M) 
NRTTE  t:6,13> 
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15     f  URMhT  vlHl. -■!■:«. r.'i.HPLUT  DK  LINITIHG  DIAMEltR  VS  IMTRUSICIN.  VQLUMt:  Q 

LT  SkhMPLE  '-.-  ^ 
Ci^LL  PLDr4  r.S.  IHTRU> 
UVITE  fb.aO) 
URITE'::iii.4> 

i\(  1  >=.; .r_4.  :>*4S4. *i;ns(s:. 'se'fie:'*.  i45>- cpr «.  oi'-'.:;4> 

Pl<IHTj-/'F'Ur-h:    iUHMLlCP    Hi    FIlLIfiii    PRE    ::URt    T:  »! .  I.i <  1  > 

c  THE  vLniniiirii,  rniipuo-   thi    M(Mt).    .(-Lanti  t-inHmr  nr-ciui    ii:e  dpiuih. 

r-  w.'.h  Vi ip I nt II  r .    HniM  iiiiii  i-i  ii'h  r  in  ciipr    jiint-uTi  p  .inn  i  iih  di    pope 

c  dihI1!.iep.     this  .  i-i.vinn  i-ihy  tt.  lh  i   duI'   unnEVLK:,   hie  huheer  400 

U  MU;T    be    IH    COl  IjMHS    l     to    3   OF    the   last    card    dp    the    PPDGRhM 

sum  ^  0.0 

-:iJM£    =    0.  0 

DO   SOO    I    ^    2.N+M+1 

nn  =     ■;Da;:'  +  d(I-i  )>-'?. 
sum  ^-  sum  -t  s;p<i>*dm 

500         SUM£   =   SUM£    ■■■    SPa>«Dt1*Dt1 

rriEHi  =  sum.--paRDi: 

liflEA-i-   =    SUri£    ■-    PDPDS 

'•.•'■i|p    =    riMEfi£    -    ririEAl»-DMErtl 

HFIEhI    =   SUMl    .'   l:USP(M+M+1> 

hi'IEh?  =  SUITS  •■'  cu:P':m+m+i> 

AVAR    =    At1EA£    -    AriEAl^AMEAl 

p:?ItiT.»=MtAN    PORE    SIZE    IS?-.  imEAl  j -'«  sSECDHD   tlOMEMr    ABOUT    THE    DRIQIM 
GIS^'DMEA£ 

PR  Wf,  ^VnR  I  (hHC  E    IS    #  >  VAR 

PkIHT./hTiJUSIEIi    ICaH    pore    size    IS.*.AriEAlj/.='HD.IU3TED    ^ECQMD  fIDMEMT 
H    AEQUT    the    origin    IS'^AMEAci 

PRIHTj^hDJUIIED    VARIAriCE    IS:»i.AVAR 

sum  --  0.0 

S;UMr:  =    0.0 

DG  >-oo  I  -  £.ri+n+i 

CM    =    <  ALDU 1  0  ■  0  ( I  >  :^  -AL  OG 1  0  •:  D  (  I  - 1  )  >  V2 . 

:.UM1    ^    sum    +    SPa>*IiM 
e.00         SUM£   =    SUM£    <-    SP<  r)*riM*DM 

Li'lEAl    =   SUML-'PORDS 

DI'1Ea£   ^   SUM,;     -   PORDS 

VhR    =    DMEAr:    --    riMEHl»Tit1EAl 

AMEAl    =^    sum     -    CUSPCH+n+l) 

AriEA£    =    SUM£    --    CMSPCN+M+l  > 

AVAR    =    HMEh£    -    AMEAl*AnEAl 

PRIHT.*MEPN   LOG    PGRE    SIZE    IS^e.  TjMEAl , /,  ;*SECDHD   MOriEMT    ABOUT    THE    DPI 
AG  IN    ISrf-riMfAS 

Ffi I rr r ,  sVhR: I hHC E    IS    5^ »  VAR 

p:;-Itir,7;rtIiJllSlf.Ii    MLHt!    LOG    PDRF    SIZE    IS^i  AHEAl » -'j  ^ADJUSTED    SECOND   nO 
PM!SHT    ABOUT    THE    OPIGItl    IS*i.AMEA£ 
400         P'<inTj*AJJJUS!l-'i:i    VhRIAHCE    IS*?AVAR 

S I  DP 

END 
C     DATA  rOLLDUS  THE  7'.--8'-9  CARD 
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c 

c 
c 
c 


u 

c 
c 
c 
c 


THIS  (--RDGRfnM  IS  ESSENTIALLY  THE  SAME:  AS  THE  PRECEEDING  PROGRAM. 

EXCEPT  IT  DHLY  PLOTS  A  SET  OF    PDRL  SIZE  CURVES  USIMG  THE  GOULD  OR 

C:^LCai1P  PLDTTLR.   FGR  E.XPLAMhTIDM  DF  THESE  PLOTTERS.  SEE  PURDUE 

DaCUMEMT  J5  CALCDMP 

pvDGRiVI  PDRE  <:  I  hi  UT  >  DUl  PUT .  TAPt:5=  I HPUT .  TAPE6=0UTPnT .  PLUT  ) 

INTEGER  TITL 

INTEGER  TITLE 

III  MEMS  I  OH  S<;5rj),pC50).D(10u)'CU<:iOO)»SP(50>fCUSP<100> 

CALL  PLOTS 

A  =  0.0 

READ.  HSET 

DO  500  K^l.MSET 

THE  PURPOSE  OF  THIS 

Di^  CURVES.   HSET  IS 

15  AN  iriTEGER  VALUE 

B  =  0.0 

CALL  FACTOR  f£.53S) 

CALCOMP  USE  CALL  FACTOR  (1 

CALL  SYMEDL^At-l.SS,  .5..7.'6 
G0.0>33> 

CALL  SYMEOLCA+l .  .  0.  75.  0.  7.--6.  . 
F^.iH.Ol       0.1        1.0 

T4E  ABOVE  PLOTS  AHD  LABELS  THE 


DO  LOOP  IS  TD 
THE  HUMEER  DF 
THAT  IS  TYPED 


PLOT  MOPE  THAM  ONE  SET 

SETS  DF  CURVES  TD  EE  PLOTTED 

ON  THE  FIRST  DATA  CARD. 


AND 


.6> 


33HLIMITIMG  PORE  DIAMETER  IM  MICROMS* 


IOC 


1000. 0.0. 0.3£) 


10.0 
X  AXIS 

CALL  SYMBOL  k A+£ . 0 . 1 . 0 . 0 . 1 . 3 . 0 . 0 . - 1 ) 
CALL  SYMBOL  ^^  ( :+3 .0.1.0.0.1.3.0.0.-1) 
C.^LL  SYMBOL  <  A+4 .  0 .  1 .  0 .  0 .  1 .  3 .  0 .  0 .  - 1 ) 
C'^LL  S'-i'MBOL  ';h+5.  0.  1.  0.  0.  1  .3.  0.  O.-l) 

TICK  MARKS  ALDHG  THE  X  AXIS 

.  1?HCUMULATIVE  POROSITY. 90. 0. 19) 

0.3H. 00. 0. 0.3) 

0.3) 

0.3) 

0.3) 


THE  ABOVE  PLACE": 
CALL  SYMBOL  'A+O 

SYMBOL 

SYMBOL 

SYMBOL 

SYMBOL 


5. 1.5- 
f.'.  .1.0. 
6.S.  0. 

(A+O. 6. 3. 0. 

\  A+O.  6.3.9. 


■vA+U 
<  A+O 


CALL 

C.^LL 

C^LL 

C.hLL 

THE  ABOVE  PLOTS  AMD 

C:-iLL  SYMBOL  <A+1 

CALL  SYMBOL  (A+1 

CALL  SYMBOL  'AH 

THE  ABOVE  PLACES 

CilLL  SYMBOL  ■' h+4 

THE  ABOVE  WRITES 


7/6. 
7.-6. 

7-6.  0.3H.  15-  0 
7-6.  0.  3H. 30. 0 
7-6.  0.3H.45.0 
LABELS  THE  Y  AXIS 
0.£. 0. 0. 1.3. 0. O.-l) 
0>3. 0. 0. 1.3. 0. O.-l) 
0.4. 0. 0. 1. 3. 0. O.-l ) 
TICK  M(^P^  S  ALDHG  THE 


CALL 
C^LL 
C'-ILL 
CALL 
CALL 
THE 


PLOT 
FLO  r 
PLOT 
PLDT 
PLOT 
ABOVE 


';a+i, 
<:a+6, 

<  A+6 , 
^A+1, 
CA+l, 


I  FGEND 
0.4.  0.3:.' 


Or.6HLEGEhD. U 
DU    THE  PLOT 


AXIS 
6) 


OF  THE  GRAPH 


PLOrTED  FOR  EACH  SET  OF 
TYPED  ON  THE  FIRST  DATA 


THE  FIRST  DATA  CARD  FDR 


0. 4. 0.2) 
0.  1.0.2) 
0.  1.0.2) 
0.4. 0.2) 

capes  drau  the  outline 
R'sad.  meat  a 
r-itiata  is  the  numbcr  df  curves  that  uill  ee 

CURVES.   MDATh  is  AN  ItlTEnEP  VALUE.  AND  IS 
CARD  FOR  THE  SET  OF  CURVES. 
DO  40  0  J=1.NDATh 

DETEPMiNATIOri  OF  PURE  SIZE  DISTRIBUTION 
PlaD.DhTE.SN 

D;i|TE  AMD  SAMPLE  NUMBER  CSN)  ARE  TYPED  DN 
EACH  CURVE. 

P'-rAD.  US.  DS.  DM.  VP,  USP.  USPM.  PE.  PF.  SR 

U:=UEIGHT  DF  SAMPLE.  DS^SPLCIFIC  GRAVITY  OF  SOLIDS.  DM=DENSITY  OF 
M::PCURY.  VP=VDLU;'1E  DF  PENETROMETER.  WSP=UEIGHT  DF  SAMPLE  AND  PEN- 
ETROMETER. USPM=IJEIGHT  DF  SAMPLE  .  PENETROMETER  AND  MFRCU9Y. 
P;£  =EVACUhTING  pressure.  PF=FILLING  PRESSURE. in  MM  DF  MERCURY. 
SR  =STEN  reading  at  filling  pressure.   THE  ABOVE  VALUES  ARE  ALL  DN 
THE  SAME  DATA  CARD 
VS=WS--DS 
UM=USPM  -  USP 
VM=UM.'-DM 

VSA=VP  -VM:*<  1 .  +PE'-PF)-SR 
VVD=VSA-VS 
VOIDP=VVD-'VS 
PORDS=VVD.'VSA 
PV=PE*VM 
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C  n  =  U:TEbEf-  VALUe  FOR  THE  NUMBEP  OF    LOU  PPECVURE  READINGS  MADE»  HOT 

C  ClUNTING  tHE  FlLl.lriG  PHEvSURE  PEmDIMG.   n=IMTLi3ER  VhLUE  FDR  THE 

C  NOMBER  Of    HIGH  PRESSURE  REhDIHGS  hflDE^  NOT  roUHTIHG  1 1  IE  IHITIftL 

C  PDRDSr  METER  PEHCiIHG.   '•  I  =  Ih  IPUSIDN  REhHING  hT  TtnTIi-iTinH  PRESSURE 

c  tH  (1L.     rnr   three  vhluk..-  hre  hll  ori  ofn    dmTh  chrd. 

C     LGU  PRESSURE  CHLCHLATIUNS 
CUO^  =  0. 
CJiPfl.'  =  0. 
DO  lOU  I=£'M+1 

p-.Hii  <  p  <:  I ;?  n  s  <  I  ■' 

C     P'l'  AND  S't::'  ARE  THE  PRESSURE  AND  INTRUSION  VALUES  FOR  EACH 

C     R£hDING>  PC  I)  IS  GIVEN  l\\    MM  DF  nERCURY  PDF  THE  LGU  PRESSURE 

C      riTRUSIDN  AND  PSI  FDR  the  HIGH  PRESSURE.   S',  I>  IS  IN  NL.   THE  TUO 

C      VALUES  ARE  TYPED  DN  DNE  DATA  CARDj  UITH  A  CARD  FDR  EACH  READING. 

FS=P' !'♦. 01934 

Vh=PV^-P':I; 

UVA=7AI-VA 

CSR=Sa)-SR 

SN=CSR-CVA 

IF  (SN.L7.fl.  ■)  GD  TD  10 

SG  =SN--US 

GSl  TD  11 

10  SG=C. 
SH  -0. 

11  cufn  =  cua-n  +  sg 

Si-'<n  =  SN'VSA 

r.jspo:'  =  cusp<i-i:'  +SP':i) 

D<  I, ■>  =  ':■ -4.  1*494.  *CDS'C£.5656>«.145)/'PS 
VhI=VA 

100   Si^-=:S':  i; 

C      HIGH  PRESSURE  CALCULSTIOtlS 

h:i  =  ii.*.  ooii»vM  -iTaiOO. 

DQ    £00    I-tU£.l1fN+l 

R.-AD'PO  ■'iSU: 

RU  '=P(I)  +  11. 

HC=R<  I  V«.  001 1*VN.-  15000. 

CMC=HC-HCI 

CSR=-;cI>  SI 

SN=CNR-CHr 

'■'-<■  1.'    =    SN -VSA 

CLispa>  =  cuspci-n  +  sp<i) 

SGH=SN/IJS 

C'jf:o=cua-i)+3GH 

rir  !,  =  <,:  _4.  :;.*434.*CnS(£.5G56;*.  145).'P(I) 

h:i=hc 
£Oij      si=s<i:) 

R'AD.    NLjN.H 
C  N.    IS    AN    iriTEbER    VALUE    UMrCH    DESIGNATES    IJHAT    SYMBOL    WILL    BE 

C  P'.DTTED    DN    THE    CURVE.       W    IS    THE    IJATER    CONTENT.       H    IS    THE    FRDST 

C  HEAVE. 

DJ  50IJ  I=£.N+i'IH 

Da- 1>  =ALnGiu';Da;o 

C      THE  PURPOSE  OF  l"HIS  DO  LDOP  IS  TD  CONVERT  THE  PORE  SIZE  DIAMETER 
C      TQ  THE  LOG  OF  THE  DIAMETER 
300    C.JSP(I-1>  =  CUSP  <  I) 

Dai+N+1;  --3.  0  -A 
C      DCM+N+rJ  IS  THF  X  hXIS  VALUE  RELATIVE  TO  THE  ORDINATE  DF  THE  PLOT 

DO^+N^£:'  =1.0 
C      D<M+H+£>  IS  THE  X  INCRFNENT  FDR  EACH  INCH  OF  PLOT 

CUSP<N+N+i:'  =-0.  15 
C      i.USP<:M+N+i::'  is  the  Y  axis  value  RELATIVE  TD  THE  DRDIMATE  OF  THE 
C      PlDT. 

CUSP<M-i-N+£:)  =0.  15 

c    CJSP<M+r^+£:■  is  the  y  increment  fdr  each  inch  dp  plot 

call  symbol  <:.  A+3 .  y£  j  3 .  7-E  >  .  0? ;  NL  j  0 .  0  n  -  n 

CALL  Si'MEDL  <A+4.  0J3.7-E.  0.  07.  15HUHTER  CONTENT  =>0.0>15) 

CALL  SYMBOL  <ft+4. O5 3. g-B) 0. 07. 19HHEHVE  AT  24  HOURS  => 0. P- 19) 

CALL  NUMBER  f A+4.  96. 3. 7-B. 0. 07. Us  0.  0.  1  > 
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C-^LL  NUMBF.K  C  A+S.  £m,  3.  t-p,  n.  07,  H.  0.  0.  1  > 
C  THE  hBQVE  lists  THE  WATER  CDHTENT  AND  FRDST  HEAVE  FDR  THE  PORE 

C      SIZE  CURVE 

criLL  Lite  'D-curpwi+n.  1,  i-riL-' 

C  THE    ABnVL    PLDT :    THE    AlTUmL    PDRL    SIZE    CURS'E 

C  I:    r    f[\F    IHE    Pl.HiFMEHr    HF    THE    FRlIil    MEhVE    AMP   UATFR   r.nrUEMT    VALUES 

■1 U  ij  }■■    -     P     ^     I.I .  '■ 

^•Cli.l  A    -    A    +    g.  0 

C      H  Iv  M-en  FDR  THE  LHCATIDH  DF  THE  POPE  SIZE  PLOT 

CALL  PLDT  '  0 .  ij  <  0 .  0 1  99'D  ) 

CTUP 

END 


Ill 


Appendix  B 

As  nentioned  in  the  body  of  the  report,  critical  region  drying  was 
tried,  but  was  unsuccessful  for  the  silt-kaolin  test  soils.   The  purpose 
of  Appendix  B  is  to  explain  the  setup  of  apparatus  and  testing  procedure 
that  v;as  used.   Explanation  of  the  soil  structural  damage  that  was 
encountered  with  this  procedure  will  be  given  in  a  future  report  by 
Ignacio  Garcia-Bengochea. 

Assembly  of  Apparatus 

The  apparatus  for  critical  region  drying  was  the  same  as  that  used 
by  Bhasin  (1975),  however  it  had  been  disassembled  and  moved.   Due  to 
reassembly,  some  leaks  were  encountered  with  the  stainless  steel  fittings. 
In  some  cases  they  could  be  eliminated  by  using  a  vise  to  hold  the 
assembly  and  then  applying  a  large  tightening  force  with  a  12  in. 
Cresent  wrench.   The  stainless  steel  fittings  are  quite  strong,  and  thus 
it  is  almost  impossible  to  overtighten.   If  this  did  not  work,  the 
fitting  was  disassembled,  the  ferrule  was  cut  off  the  tube,  a  new  ferrule 
was  inserted,  and  the  fitting  was  tightened.   Once  all  leaks  had  been 
corrected,  the  system  was  able  to  maintain  a  pressure  of  4000  psi,  with 
pressure  changes  due  solely  to  temperature  fluctuations. 

Along  with  reassembly  of  the  critical  region  device,  a  shield  was 
constructed  to  protect  against  any  leaking  superheated  steam.   The  shield 
was  constructed  of  0.125  in.  mild  steel  plates  bolted  to  a  angle  iron 
framework.   The  shield  was  designed  so  that  a  cover  plate  could  be 
removed  to  insert  the  critical  region  drying  vessel  into  the  furnace. 
The  cover  plate  was  left  off  the  shield  during  the  early  part  of  the 
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critical  region  run  in  order  to  check  for  leaks.   Once  the  pressure  had 
risen  to  about  500  psi,  the  cover  plate  was  bolted  to  the  rest  of  the 
shield. 

Procedure 
The  procedure  used  during  a  critical  region  run  was  essentially  the 
samj  as  that  used  by  Bhasin  (1975).   However,  the  follov/ing  points  are 
listed  in  order  to  clarify  the  method. 

1)  The  0.5  in.  diameter  samples  were  trimmed  using  a 
push  tube  similar  to  that  used  by  Bhasin,  however, 
it  was  constructed  so  that  the  sharpened  tube  could 
be  replaced.   Also,  the  sampler  was  designed  so  that 
the  sample  container  fitted  into  it,  and  thus  the 
sample  could  be  Intruded  directly  into  the  container. 
Once  the  sample  was  trimmed,  the  container  and  sample 
were  removed  from  the  sampler.   The  sample  ends  were 
then  trimmed  flush  to  the  0.5  in.  high  sample  container 
using  a  razor  blade. 

2)  The  carborundum  porous  stones  that  were  placed  on  each 
end  of  the  sample  container  are  not  commercially  available 
for  the  small  diameter  needed.  Therefore,  they  were 
broken  from  large  diameter  0.25  in.  thick  porous  stones 
with  a  cold  chisel  and  pilars  to  a  size  to  cover  the 
sample  container.   It  was  not  necessary  that  the  stones 

be  round,  only  that  they  fit  between  the  stainless  steel 
'    end  plates  and  the  sample  container. 
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3)  Following  assembly  of  the  sample  containers,  they 
were  soaked  in  water  till  the  critical  region  run 
was  initiated.   Since  the  samples  were  trimmed 
from  soil  that  was  saturated,  it  was  not  necessary 
to  saturate  the  small  samples  as  Bhasin  did. 

4)  The  samples  were  placed  in  the  critical  region  vessel, 
and  water  was  added  while  the  vessel  sat  in  the  steel 
holder  attached  to  the  wall.  The  vessel  top  was 
lubricated  using  Never  Seez,  and  then  screwed  onto 
the  body  of  the  vessel.   The  top  was  secured  by 
tightening  the  six  screws  in  a  staggered  order.   The 
screws  were  tightened  to  70  ft  lbs  with  a  torque  v/rench. 

5)  The  vessel  was  placed  inside  the  furnace  and  the 
thermocouple  wire  was  attached.   Two  asbestos  cover 
plates  were  placed  on  top  of  the  furnace  for 
insulation  and  the  pressure  line  was  attached  to  the 
main  pressure  system. 

6)  Following  connection  of  the  vessel,  the  pressure  was 
raised  to  200  psi  with  the  hydraulic  jack.   To  use 
the  jack,  the  valve  from  the  jack  to  the  system 

was  closed,  after  which  the  valve  to  the  plastic 
vrater  vessel  was  opened.   To  allow  water  to  enter 
the  jack,  it  was  jacked  outward,  following  which 
the  valve  to  the  water  vessel  was  closed.   The 
pressure  was  then  raised  to  the  system  pressure,  and 
the  valve  leading  to  the  system  was  opened. 
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7)  The  temperature  and  pressure  during  the  critical  region 
run  were  raised  using  the  time  relationship  shown  in 
Figure  39.  The  heater  temperature  was  controlled  using 
a  Research  Incorporated  Thermae  Temperature  Controller 
on  the  "set  point"  mode.  The  controller  was  calibrated 
using  an  ice  bath,  boiling  water,  and  the  melting  point 
of  lead. 

8)  To  raise  the  temperature,  the  temperature  selector  on 
the  Thermae  Controller  was  set  about  15  F  above  the 
current  temperature  and  the  rate  of  heating  was 
adjusted  using  the  proportional  band  control  so  that 
about  200  to  800  watts  were  applied  to  the  heater. 

The  heater  wattage  was  measured  using  a  watt-amp  meter 
connected  to  the  temperature  controller. 

9)  For  the  critical  region  run  up  to  about  2000  psi,  the 
pressure  was  controlled  by  using  a  nitrogen  cylinder 
and  regulator.   Above  2000  psi,  the  nitrogen  tank  was 
shut  off  from  the  system  and  the  pressure  was  controlled 
by  the  hydraulic  jack.   Because  of  the  effects  of 
rising  temperature  on  pressure,  the  pressure  usually  had 
to  be  backed  off.   The  jack  could  be  emptied  or  filled 
during  the  run  using  the  earlier  mentioned  procedure. 

10)  Once  the  critical  region  run  was  completed,  and  the 
pressure  was  back  to  zero,  the  vessel  pressure  tube  was 
disconnected  from  the  rest  of  the  system  at  the  first  fitting. 

11)  Following  disconnection,  the  vessel  was  cooled  overnight  and 
then  disassembled. 
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Appendix  C 


Purdue  University  Negative  Numbers  for  Photographs 


Figure  Number 

2 

3 

5 

6 

7 

8 

9 
10 
11 
12 
13 
14 
22 


23 


Negative  Number 


Dry  of  Optimum 

Optimum 

Wet  of  Optimum 

Dry  of  Optimum 

Optimum 

Wet  of  Optimum 


75529 

24 

75529 

34 

75529 

32 

75609 

12 

75529 

30 

75609 

16 

75609 

19 

75529 

18 

75529 

25 

75529 

28 

75609 

14 

75609 

5 

75529 

1 

75529 

4 

75529 

3 

75529 

9 

75529 

10 

75529 

13 
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tr 

LU 
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